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SECTION  I 
INTRODUCTION 


This  report  provides  user  instructions  for  the  1986  version  of  the  EPIC-2  code. 

The  first  documented  version  of  the  EPIC-2  code  was  developed  for  the  Ballistic  Re¬ 
search  Laboratory  (BRL)  in  1977  (Reference  1).  The  second  documented  version  was 
developed  for  the  Air  Force  Armament  Laboratory  (AFATL)  in  1983  (Reference  2). 
The  1983  version  added  expanded  preprocessing  and  postprocessing  capabilities,  a 
heat  conduction  option,  improved  material  models,  a  material  data  library,  im¬ 
proved  sliding  interface  capabilities,  and  a  more  accurate  mass-lumping  for 
mulation. 

The  1986  version,  described  in  this  report,  incorporates  a  NABOR  option  for  vari 
able  nodal  connectivity.  The  essential  nature  of  this  option  is  that  each  node  is  af 
fected  only  by  its  nearest  neighbor  nodes.  As  the  nodes  move  closer  than  their 
equilibrium  distance,  they  generate  compressive,  repulsive  forces.  Conversely,  when 
they  move  apart,  they  generate  tensile,  attractive  forces.  Material  strength  effects 
are  also  included.  The  key  to  this  approach  is  that  it  is  possible  to  have  variable  con¬ 
nectivity;  i.e.,  a  node  can  take  on  new  nearest  neighbors,  thus  allowing  all  forms  of 
distortion.  Based  on  the  concept  of  nearest  neighbors,  this  option  has  been  desig 
nated  the  NABOR  option.  A  complete  description  of  this  approach  will  be  presented 
in  a  forthcoming  final  report  for  this  contract. 

Other  new  features  of  the  1986  version  are  capabilities  to  include  crushable  materi 
nls  •'such  as  concrete),  one  dimensional  geometry,  plane  stress  geometry,  two- 
dimensional  quad  elements  (in  addition  to  the  standard  triangular  elements),  erod 
ing  interfaces  for  projectile  penetration  into  thick  targets,  expanded  postprocessing, 
and  a  target  add  option. 

Although  quad  elements  are  being  introduced  as  an  option,  one  of  the  distinguishing 
features  of  the  EPIC-2  code  is  that  it  is  primarily  based  on  a  triangular  element  for 
mulation.  This  allows  problems  with  severe  distortions  to  be  run  without  rezoning. 
This  is  a  very  desirable  feature  for  many  intense  impulsive  loading  problems  involv 
ing  high  velocity  impact  or  explosive  detonation. 


SECTION  II 
USER  INSTRUCTIONS 


This  section  provides  user  instructions  for  the  EPIC-2  code.  The  code  consists  of  a 
Preprocessor,  Main  Routine  and  Postprocessor  for  state  and  time  plots.  The  for¬ 
mulation  is  not  provided  here;  however,  most  of  the  equations  are  identical  to  those 
of  the  original  version,  (Reference  1).  Some  comments  concerning  the  formulation  of 
the  current  version  of  the  code  follow: 

•  The  crossed  triangle  arrangement  of  triangular  elements  (four  triangles  in  a 
quadrilateral)  has  been  shown  to  give  better  accuracy  and  should  be  used  in 
most  instances  (Reference  3). 

•  The  sliding  interface  formulation  is  similar  to  that  used  in  the  three- 
dimensional  EPIC-3  code  (References  4  and  5). 

•  The  eroding  interface  option  is  described  in  Reference  6. 

•  The  heat  conduction  option  is  described  in  Reference  7. 

•  The  plane  stress  option  is  described  in  Reference  8. 

•  The  NABOR  and  quad  element  formulations  will  appear  in  the  final  report  for 
this  contract.  The  NABOR  formulation  for  plane  strain  is  given  in  Reference  9. 

•  A  characteristic  of  triangular  elements  is  that  they  sometimes  produce  a  pres 
sure  field  that  oscillates  spatially.  This  can  cause  inaccuracies  for  cases  where 
the  material  strength  or  fracture  is  dependent  on  the  pressure  in  the  element. 
This  is  overcome  by  computing  a  nodal  pressure  which  is  the  average  of  the  ele¬ 
ment  pressures  for  all  elements  which  contain  a  specified  node.  The  nodal  pres¬ 
sures  do  not  oscillate  spatially  but  rather  form  a  smooth  pressure  field.  A 
smoothed  eler.  ent  pressure  is  then  defined  as  the  average  of  the  three  as¬ 
sociated  nodal  pressures.  This  smoothed  pressure  is  used  for  the  strength  and 
fracture  models  which  are  pressure  dependent  (Reference  3). 

A  description  of  input  data  for  the  EPIC-2  code  is  given  in  Figures  1  through  15. 
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INPUT  DATA  FOR  THE  PREPROCESSOR 


The  function  of  the  Preprocessor  is  to  define  the  initial  geometry  and  velocity  condi 
tions.  The  descriptions  which  follow  are  for  the  data  in  Figure  1.  Consistent  unit: 
must  he  used  and  the  unit  of  time  must  he  seconds. 


Description  Card  (2I5.A60)  — 

CASE  =  Case  number  for  run  identification. 


TYPE 


1  specifies  a  Preprocessor  run  only. 

2  specifies  a  Preprocessor  and  Main  Routine  Run. 


DESCRIPTION  =  Description  provided  by  user. 
OF  PROBLEM 


Miscellaneous  Card  (1615)  — 

II  specifies  axisymmetric  geometry 
2  specifies  plane  strain  geometry 
3  specifies  one-dimensional  geometry  (cartesian) 

4  specifies  plane  stress  geometry 

10  will  not  print  individual  data  for  each  node  and  element 

1  will  print  individual  data 

10  will  not  write  Preprocessor  data  on  restart  tape 

1  will  write  Preprocessor  data  on  restart  tape 

NMAS  Number  of  concentrated  masses  to  be  input  separately. 

NRIG  Number  of  systems  of  nodes  which  move  as  rigid  system. 

NCHNK  -  Number  of  groups  or  elements  for  which  subsystem  data  are  re 

quested. 


NSLID 


Number  of  sliding  surfaces. 


T? 


1  gives  a  rigid  frictionless  surface  on  the  positive  side  of  the 
plane  described  by  z  =  0.  If  the  equations  of  motion  cause  a  node 
to  have  a  negative  z  coordinate,  the  z  coordinate  and  velocity  are 
set  to  zero.  If  IZR  =  0  this  option  is  not  used. 


SPLIT 


0  will  perform  the  sliding  surface  computations  after  the  up¬ 
dated  velocities  and  displacements  are  determined  from  the  usu¬ 
al  equations  of  motion.  Contact  ^established  as  long  as  the 
slave  node  interferes  with  the  master  surface  before  the  veloc¬ 
ities  and  displacements  are  adjusted  in  the  sliding  surface  rou 
tines.  This  option  is  the  most  reliable  and  should  be  used  for 
complicated  sliding  surfaces  which  include  double  pass  options 
and  intersecting  sliding  surfaces.  It  must  be  used  for  the  eroding 
interface  option.  For  relatively  low  velocity  impact  problems, 
where  there  is  limited  deformation,  this  option  can  introduce 
significant  errors  in  the  form  of  excessive  deformation  and  in 
ternal  energy. 

1  will  perform  the  sliding  surface  computations  after  the  up 
dated  velocities  are  determined  but  before  the  updated  displace¬ 
ments  are  determined.  This  should  be  used  if  the  sliding  surfac¬ 
es  are  relatively  simple.  Contact  is  first  established  when  the 
slave  node  interferes  with  the  master  surface.  Thereafter,  a 
slave  node  is  considered  to  be  in  contact  until  the  preadjusted 
normal  velocities  between  the  two  surfaces  are  separating  rath 
er  than  closing.  This  approach  minimizes  the  distance  the  slave 
node  is  moved  to  place  it  on  the  master  surface  and  is  therefore 
more  accurate.  If  there  are  no  sliding  surfaces  either  option  can 
be  used  (SPLIT  =  0  or  SPLIT  =  1). 


NTOP 


Node  number  of  the  inner  radii  node  at  the  top  free  end  of  a 
projectile.  These  data,  together  with  NBOT  and  NRING,  are 
used  if  it  is  desired  to  obtain  internal  loads  in  a  slender 
projectile.  Use  primarily  for  axisymmetric  geometry  (GEOM 
1)  and  one-dimensional  geometry  (GEOM  =  3).  It  can  also  be 
used  for  normal  (but  not  oblique)  impact  along  the  Z  axis  with 
plane  strain  or  plane  stress  geometry  (GEOM  =  2  or  4). 


NBOT 


Node  number  of  the  inner  radii  node  at  the  lower  end  of  the 
projectile.  Use  only  if  internal  loads  are  desired. 


NRING 


Number  of  rings  of  nodes  in  the  projectile.  Includes  the  secon 
dary  nodes  for  crossed  triangle  geometry.  The  nodal  geometry 
must  be  consistent  with  that  of  the  rod  shape  geometry  shown  in 
Figure  10.  Use  only  if  internal  loads  are  desired. 
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AVEP 


CDUCT 


UNIT 


NRST 


0  will  compute  individual  pressures  for  each  element.  This  is  the 
primary  option  and  should  be  used  with  triangular  elements  for 
problems  involving  severe  distortions  or  explosive  material'.  It 
is  also  used  for  quad  elements  and  one  dimensional  geometry 

2  will  compute  an  average  pressure  for  each  group  of  two  trian¬ 
gular  elements.  The  elements  are  grouped  in  pairs  based  on  the 
order  in  which  they  are  input.  Damage  can  be  computed,  but 
fracture  is  not  allowed.  Only  solid  materials  (no  liquids,  ex¬ 
plosives  or  crushable  materials)  and  triangular  elements  (no 
quads)  can  be  used. 

4  will  compute  an  average  pressure  for  each  group  of  four  trian¬ 
gular  elements  (crossed  triangle  geometry)  similar  to  that  de¬ 
scribed  above.  It  should  be  used  only  for  cases  in  which  the  pres¬ 
sure  is  not  significantly  greater  than  the  strength  of  the  materi¬ 
al. 

Note  for  AVEP  =  2  or  4:  With  the  introduction  of  quad  ele¬ 
ments,  the  primary  use  for  AVEP  =  2  or  4  is  when  the  heat  con¬ 
duction  option  is  used  (CDUCT  =1).  When  the  heat  conduction 
option  is  not  used,  quad  elements  are  recommended  instead  of 
the  AVEP  =  2  or  4  option. 

1  incorporates  the  heat  conduction  option.  If  CDUCT  =  1  then 
the  average  pressure  option  must  be  used  (AVEP  =  2  or  4).  Heat 
conduction  is  not  included  for  CDUCT  =  0. 

0  indicates  the  constants  in  the  material  library  have  English 
units  (pound/inch/second/degree  Fahrenheit). 

■  1  indicates  the  constants  in  the  material  library  are  converted  to 
Standard  International  (SI)  units. 


=  Number  of  groups  of  nodes  to  have  restraints  redefined. 


Material  Data  Cards  —  Material  data  can  be  completely  defined  by  the  user  or 
taken  from  the  material  data  library.  Specific  instructions  are  presented  later.  End 
material  data  with  a  blank  card. 


Projectile  Scale/Shift/Rotate  Card  (5F10.0)  — 


RSCALE  -  Factor  by  which  the  r  coordinates  of  all  projectile  nodes  are  mul 

tiplied.  Applied  after  the  coordinate  shifts  (RSHIFT,  ZSHIFT) 
and  before  the  rotation  (ROTATE)  described  later.  Must  use 
RSCALE  -  1.0  if  NABOR  nodes  are  used. 

ZSOALE  -  Factor  by  which  the  z  coordinates  are  multiplied.  Must  use 

ZSCALE  1.0  if  NABOR  nodes  are  used 

ROTATE  Rotation  about  r  z  =  0  in  the  r-z  plane  of  all  projectile  nodes 

(degrees)  (positive  clockwise).  Applied  after  the  coordinate  shifts 
(RSHIFT,  ZSHIFT)  and  the  scale  factors  (XSCALE,  ZSCALE i. 

RSHIFT  Increment  added  to  the  r  coordinates  of  all  projectile  nodes 

(length).  Applied  before  the  scale  factors  (RSCALE,  ZSCALE). 

ZSHIFT  -  Increment  added  to  the  z  coordinates  (length). 

Node  Data  Cards  for  Projectile  —  These  cards  are  required  to  define  the  projectile 
nodes.  If  a  node  is  at  the  interface  of  the  projectile  and  the  target  and  contains  mass 
from  both  the  projectile  and  the  target,  it  must  be  included  with  the  projectile  nodes. 
The  node  numbers  must  not  exceed  the  dimension  of  the  node  arrays,  and  they  need 
not  be  numbered  consecutively  or  in  increasing  order.  Specific  instructions  for  node 
input  data  are  presented  later.  End  projectile  node  data  with  a  blank  card. 

Target  Scale/Shift/Rotate  Card  (5F10.0)  —  Same  as  Projectile  Scale/Shift/Rotate 
Card  except  it  applies  to  the  target  nodes.  Must  be  included  even  if  there  are  no 
target  nodes. 

Node  Data  Cards  for  Target  —  Similar  to  node  data  cards  for  projectile.  Specific 
instructions  are  presented  later.  End  target  node  data  with  a  blank  card.  Include 
blank  card  even  if  there  are  no  target  nodes. 

Element  Data  Cards  for  Projectile  —  These  cards  are  required  to  define  the 
projectile  elements.  The  element  numbers  must  not  exceed  the  dimension  of  the  ele¬ 
ment  arrays,  and  they  need  not  be  numbered  consecutively  or  in  increasing  order. 
Specific  instructions  are  presented  later.  End  projectile  element  data  with  a  blank 
card. 

Element  Data  Cards  for  Target  —  Similar  to  element  data  cards  for  projectile. 
Specific  instructions  are  presented  later.  End  target  element  data  with  a  blank  card 
Include  blank  card  even  if  there  are  no  target  elements. 


Concentrated  Mass  Cards  (15,  5X,  F10.0)  —  There  are  NMAS  (defined  in  Mis 
cellaneous  Card)  cards  entered  for  the  concentrated  masses.  These  cards  are  omitted 
when  NMAS  =■  0.  Each  card  contains  data  for  one  mass.  The  program  is  currently 
limited  to  50  concentrated  masses.  Concentrated  masses  cannot  be  added  to  NABOR 
nodes. 


N  -  Node  number  to  which  the  concentrated  mass  is  added 

MASS  (N)  =  Concentrated  mass  added  to  node  N 

Rigid  Body  Identification  Cards  (215)  —  Each  system  of  rigid  body  nodes  contains 
one  Rigid  Body  Identification  Card  and  additional  cards  to  specify  the  nodes.  The 
program  is  currently  dimensioned  for  a  maximum  of  10  rigid  body  systems  with  a 
maximum  of  50  nodes  per  system.  If  there  are  no  rigid  body  systems  (NRIG  =  0  in 
Miscellaneous  Card)  this  group  of  cards  is  omitted.  If  there  is  more  than  one  system 
of  rigid  body  nodes,  all  cards  for  the  first  system  are  entered  before  the  Rigid  Body 
Identification  Card  for  the  next  system  is  entered.  Rigid  body  nodes  must  not  con¬ 
tain  any  slave  nodes  on  sliding  interfaces  or  have  nodes  restrained  in  the  z  direction. 
For  axisymmetry  geometry  (GEOM  =  1),  all  rigid  body  nodes  are  restrained  in  the  r 
direction.  For  plane  strain  or  plane  stress  geometry  (GEOM  =  2  or  4),  if  any  rigid 
body  node  is  restrained  in  the  r  direction,  then  all  are  restrained  in  the  r  direction. 

NRN  =  Number  of  rigid  body  nodes  in  the  system 

NRG  -  Number  of  groups  of  rigid  body  nodes  to  be  read.  If  NRG  =  0  the 

nodes  are  read  individually. 

Individual  Rigid  Body  Nodes  Cards  (1615)  — 

R1  RN  --  Individual  nodes  in  rigid  body  system 

NRG  Grouped  Rigid  Body  Nodes  Cards  (315)  — 

RIG  First  node  in  the  group  of  rigid  body  nodes 

RNG  Last  node  in  the  group  of  rigid  body  nodes 

INC  Increment  between  nodes  in  the  group  of  rigid  body  nodes 


Chunk  Data  Identification  Cards  (215)  —  Each  subsystem  of  element  chunks  for 
which  output  is  desired  requires  a  Chunk  Data  Identification  Card  and  additional 
cards  to  specify  the  specific  elements.  The  program  is  currently  dimensioned  for  a 
maximum  of  10  subsystems  of  chunks  with  a  maximum  of  50  elements  per 
subsystem.  If  there  are  no  chunk  data  to  be  obtained  (NCHNK  =  0  in  miscellaneous 
card)  this  group  of  cards  is  omitted.  Data  are  input  in  a  manner  similar  to  that  of  the 
rigid  body  systems  described  previously. 

NCE  =  Number  of  elements  in  the  subsystem 

NEC  Number  of  groups  of  elements  to  be  read 

Individual  Element  Chunk  Cards  (1615)  — 

E1...EN  =•-  Individual  elements  in  subsystem 

NEC  Grouped  Element  Chunk  Cards  (315)  — 

ElG  =  First  element  in  the  group  of  chunk  elements 

ENG  -  Last  element  in  the  group  of  chunk  elements 

INC  Increment  between  elements  in  the  group  of  chunk  elements 

Slide  line  Identification  Cards  (815,  3FI0.0)  —  Each  sliding  surface  contains  one 
Identification  Card  and  cards  (as  required)  describing  the  master  nodes  and  slave 
nodes.  The  program  is  dimensioned  for  a  maximum  of  10  sliding  surfaces  which  can 
each  contain  a  maximum  of  50  master  nodes  and  500  slave  nodes.  If  there  are  no 
sliding  surfaces  (NSLID  -  0  on  Miscellaneous  Card),  this  group  of  cards  is  omitted. 
If  there  is  more  than  one  sliding  surface,  all  data  for  the  first  sliding  surface  are 
entered  before  the  second  Sliding  Surface  Identification  Card  is  entered.  The  master 
nodes  generally  should  include  the  higher-density  material.  It  is  also  desirable  for 
the  master  surface  to  have  the  stronger  material,  have  equal  or  greater  spacing  than 
the  slave  nodes,  and  not  have  a  convex  surface  toward  the  slave  nodes. 

NMN  -  Number  of  master  nodes  on  the  sliding  surface 

NSN  Number  of  slave  nodes  on  the  sliding  surface.  If  NSN  ~  0,  then 

slave  nodes  will  Ik*  specified  by  the  region  using  the  Slave  Node 
I  .imils  ( 'ai  d. 

NMG 


NSG 


Number  of  groups  of  master  nodes  to  be  read 
Number  of  groups  of  slave  nodes  to  he  read 


m 


Number  of  velocity  iterations  for  the  slave  nodes  on  the  master 
surface  (References  4  and  5).  Errors  in  the  velocity  match  lead  to 
errors  in  the  deviator  and  shear  stresses,  but  generally  not  the 
pressure.  For  high-velocity  impact  and  explosive  detonation, 
where  the  pressures  are  much  higher  than  the  deviator  and 
shear  stresses,  a  relatively  low  value  of  IT1  =  1  or  ITl  =  2  can 
be  used.  For  lower  pressure  problems,  higher  values  should  be 
used,  ITl  =  2  to  ITl  =  5.  The  velocity  iterations,  and  the  cor 
responding  searches  on  the  master  surface,  are  performed  only 
for  those  slave  nodes  found  to  be  in  contact  during  the  first  iter¬ 
ation.  For  one-dimensional  geometry  (GEOM^3)  and  the  eroding 
interface  option  (ERODE  >  0),  use  ITl  -  1. 

1T2  -  Number  of  velocity  iterations  of  the  master  nodes  on  the  slave 

surface.  This  allows  a  double  pass  to  he  made  such  that  there  is 
no  interference  or  crossover  on  the  sliding  surface  (References  4 
and  5).  If  IT2  =  0,  there  is  no  second  pass,  and  the  slave  nodes 
can  be  input  in  any  order.  With  this  option  it  is  possible  to  desig¬ 
nate  interior  nodes  (as  well  as  surface  nodes)  as  slave  nodes. 
This  procedure  allows  elements  containing  slave  nodes  to  fail 
completely  to  simulate  an  eroding  process  (Reference  4).  For  IT2 
>  0,  a  double  pass  is  performed,  and  the  slave  nodes  must  be  in 
put  in  a  specified  order.  The  double  pass  option  (IT2  >  0)  can 
only  be  used  with  SEEK  =  1. 

II  is  the  option  used  for  two-dimensional  (GEOM  -  1,2,4)  prob¬ 
lems  with  no  NABOR  slave  nodes. 

2  is  the  option  used  when  all  slave  nodes  are  NABOR  nodes 
(NABOR  nodes  cannot  be  master  nodes). 

3  is  the  option  used  for  one-dimensional  geometry  (GEOM  =  3). 

MBOT  -  Lowest  number  node  on  the  bottom  of  the  target  plate  when  per 

forming  perforation  computations  with  the  eroding  interface  op¬ 
tion  (ERODE  >  0).  Ail  nodes  above  the  bottom  surface  of  the 
plate  must  have  lower  node  numbers  than  MBOT.  This  criterion 
is  satisfied  when  using  the  flat  plate  generator  for  nodes  and  ele 
ments. 

FRICTION  -  Coefficient  of  sliding  friction  in  the  r-z  plane.  Does  not  act  in  0 
direction  for  relative  spinning  velocities. 
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HE?'  VEL  =  Reference  velocity,  which  when  multiplied  by  the  integration 

time  increment,  gives  a  reference  distance.  Slave  nodes  are  con¬ 
sidered  to  be  associated  with  a  particular  master  surface  only 
when  they  are  within  this  reference  distance.  It  is  recommended 
that  REF  VEL  be  1.0  to  1.5  times  the  initial  relative  impact  ve¬ 
locity  or  the  detonation  velocity  of  explosives  contained  in  the 
problem. 

ERODE  =  Erosion  strain  for  the  eroding  interface  option.  This  option  is  for 

penetration/perforation  of  thick  plates  and  is  activated  when 
ERODE  >  0.  It  applies  to  triangular  elements  only  and  should 
only  be  used  when  erosion  is  the  primary  mode  of  penetration. 
The  algorithm  and  example  problems  are  presented  in  Reference 
6.  Because  total  failure  of  the  elements  must  be  achieved  by  the 
eroding  interface  algorithm,  it  is  important  that  EFAIL  (a  mate¬ 
rial  property)  be  much  greater  than  ERODE. 

An  eroding  interface  usually  consists  of  two  sets  of  slide  line 
data.  The  first  slide  line  usually  designates  the  top  surface  of  the 
plate  as  the  master  surface  and  the  potentially  eroded  nodes  in 
the  projectile  as  slave  nodes.  The  second  slide  fine  usually  desig 
nates  the  outer  surface  of  the  projectile  as  the  master  surface 
and  the  potentially  eroded  nodes  in  the  plate  as  slave  nodes. 
This  ensures  that  there  is  no  cross-over  of  material  between  the 
projectile  and  the  target.  Under  some  instances  it  may  only  be 
possible  to  use  the  first  slide  line. 

When  using  the  erosion  option  with  axisymmetric  geometry 
(GEOM  ~  1),  the  radial  velocities  of  free  nodes  (those  whose  as¬ 
sociated  elements  have  all  eroded)  are  adjusted  based  on  the  cor 
responding  axial  velocity  (RDOT  =  .5|ZDOT|).  Therefore,  the  ini¬ 
tial  target  velocity  in  the  axial  direction  should  be  zero, 
otherwise  the  radial  velocities  would  be  adjusted  incorrectly. 

Individual  Master  Node  Cards  (1615)  —  The  master  nodes  are  input  individually 
when  NMG  =  0.  The  nodes  must  be  entered  in  order  from  the  first  master  node  to 
the  last  master  node  along  the  row  of  nodes.  When  moving  from  the  first  node  to  the 
last  node,  the  slave  nodes  must  be  to  the  left  of  the  master  surface.  For  one¬ 
dimensional  geometry  (GEOM  -  3),  the  master  node  must  be  below  (lower  z  coordi 
nate)  the  associated  slave  node. 

Ml  ,MN  Master  nodes  in  proper  order  from  Ml  to  MN 
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NMG  Grouped  Master  Node  Cards  (315)  —  This  option  allows  master  nodes  to  lie 
input  in  groups.  However,  the  nodes  must  be  input  in  the  proper  order,  as  described 
Cor  the  Individual  Master  Node  Cards. 

MIG  =  First  node  in  the  group  of  master  nodes. 

MNG  -  Last  node  in  the  group  of  master  nodes. 

INC  --  Increment  between  the  nodes  in  the  group  of  master  nodes 

Individual  Slave  Node  Cards  (1615)  —  The  slave  nodes  are  input  individually 

when  NSG  =  0.  The  slave  nodes  must  be  to  the  left  of  the  master  surface  when  mov 
ing  from  the  first  master  node  to  the  last  master  node.  If  there  is  no  double  pass  (IT2 
=  0),  the  slave  nodes  can  be  input  in  any  order.  If  the  double  pass  option  is  used  (IT2 
>  0),  the  slave  nodes  must  be  input  in  the  opposite  direction  of  the  master  nodes. 
The  first  slave  node  must  be  near  the  last  master  node  and  the  last  slave  node  must 
be  near  the  first  master  node.  This  means  the  master  surface  is  to  the  left  of  the 
slave  surface  when  moving  from  the  first  slave  node  to  the  last  slave  node.  For  one 
dimensional  geometry  (GEOM  -  3),  the  slave  node  must  be  above  (higher  z  roordi- 
ante)  the  associated  master  node. 

S1...SN  =  Slave  nodes  in  proper  order  from  Si  to  SN 

NSG  Grouped  Slave  Node  Cards  (315)  —  This  option  allows  the  slave  nodes  to  be 
input  in  groups.  Restrictions  on  order  of  input  are  as  described  for  the  Individual 
Slave  Node  Cards. 

S1G  First  node  in  the  group  of  slave  nodes. 

SNG  -  Last  node  in  the  group  of  slave  nodes. 

INC  =  Increment  between  nodes  in  the  group  of  slave  nodes. 

Slave  Node  Limits  Card  —  For  NSN  =  0  (4F10.0)  —  This  option  allows  all  nodes 
within  a  specified  region  to  be  designated  as  slave  nodes.  Used  only  if  NSN  =  0  on 
Slide  Line  Identification  Card. 

RMAX  -  Maximum  r  coordinate  of  slave  node  region. 

RMIN  =  Minimum  r  coordinate  of  slave  node  region. 

ZMAX  Maximum  z  coordinate  of  slave  node  region. 

-  Minimum  z  coordinate  of  slave  node  region. 


ZMIN 


Restrained  Nodes  Identification  Cards  (315)  —  Each  set  of  restrained  nodes  con 
tains  one  Restrained  Nodes  Identification  Card  and  additional  cards  to  specify  the 
nodes.  The  program  does  not  impose  any  constraint  on  the  number  of  sets  and  each 
set  can  contain  as  many  as  the  node  arrays  can  handle.  If  there  are  no  restrained 
node  sets  (NRST  =  0  in  Miscellaneous  Card)  this  group  of  cards  is  omitted.  If  there  is 
more  than  one  set  of  restrained  nodes,  all  cards  for  the  first  set  are  entered  before 
the  Restrained  Nodes  Identification  Card  for  the  next  set  is  entered.  This  input  re¬ 
defines  the  restraints  on  the  designated  nodes  lit  does  not  simply  add  to  existing  re¬ 
straints). 

NFN  -  Number  of  nodes  in  set. 

NFG  =  Number  of  groups  of  nodes  to  be  read.  If  NFG  0  the  nodes  are 

read  individually. 

1R,  IZ,  IT  r  1  restrains  nodes  in  r,  z,  ft  directions,  respectively.  Expanded  de 

scription  given  for  Line  of  Nodes  Description  Card  in  Node 
Geometry  Subsection. 

Individual  Restrained  Nodes  C'ards  —  For  NFG  =  0(1615)  — 

FI  ,FN  Individual  nodes  to  be  restrained. 

NFG  Group  Restrained  Nodes  Cards  (315)  — 

FIG  -  First  node  in  the  group  of  nodes  to  be  restrained. 

FNG  Last  node  in  the  group  of  nodes  to  be  restrained.  (May  be  zero  if 

there  is  only  one  node  in  the  group.) 

INC  Increment  bet  ween  nodes  in  the  group  of  restrained  nodes.  (May 

be  zero  if  the  increment  is  one  or  if  no  increment  is  needed.) 

Velocity/Detonation  Card  (8F10.0)  —  This  card  describes  the  initial  velocity 
and/or  detonation  conditions.  If  there  are  interface  nodes  which  include  mass  from 
Goth  the  projectile  and  the  target,  the  velocities  of  these  nodes  are  automatically  ad 
justed  to  conserve  momenta. 

PRI)OT  Projectile  velocity  in  the  r  direction  (distance/time). 

PZDOT  Projectile  velocity  in  the  z  direction. 

IT  DOT  Projectile  velocity  in  the  0  direction  (radians/time).  Must  be  zero 

if  any  quad  elements  or  NAHOR  nodes  are  used. 


TZIX)T  Target  velocity  in  the  z  direction.  Should  be  zero  when  the 

erosion  option  is  used  (ERODE  >0)  with  axisymmetric 
geometry. 

TTDOT  =  Target  velocity  in  the  6  direction.  Must  be  zero  if  any  quad  ele 

inents  or  NABOR  nodes  are  used. 

KDET  Radial  coordinate  of  the  explosive  detonation  (distance! 

ZDET  =  Axial  coordinate  of  the  explosive  detonation 

Initial  Integration  Card  (8F10.0)  — 

DTI  =  Integration  time  increment  for  the  first  cycle.  This  must  be  less 

than  the  time  required  to  travel  across  the  minimum  dimension 
of  each  element  at  the  sound  speed  of  the  material  in  that  ele¬ 
ment. 

TBURN  =  Time  at  which  the  detonation  begins  at  RDET,  ZDET 

BETA  =  Ratio  of  compressed  time  to  actual  time.  This  is  used  to  simulate 

quasi-static  or  relatively  slow  dynamic  problems  in  which  the 
problem  must  be  run  on  a  faster  time  scale.  The  strain-rate-de 
pendent  material  models  and  heat  conduction  equations  include 
the  effect  of  BETA.  See  Reference  7  for  more  details.  Set  BETA 
=  1.0  for  normal  conditions. 

Material  Cards  for  Solids  from  Library  (415,  2F5.0)  —  Data  for  some  materials 
are  available  from  the  material  library  in  subroutine  MATLIB.  The  specific  materi¬ 
als  are  listed  as  output  from  the  Preprocessor.  Library  materials  may  be  used  direct¬ 
ly  without  being  called  by  this  card.  If  they  are  called  by  this  card,  however,  the  ma¬ 
terial  data  will  be  printed  as  part  of  the  output.  This  card  also  allows  the  user  to 
make  some  decisions  regarding  fracture  options.  The  user  should  read  the  comments 
in  subroutine  MATLIB  to  obtain  the  references  from  which  the  data  were  generated. 

MATL  =  Material  identification  number.  It  must  be  in  the  range  of  1 

through  50  and  must  correspond  to  a  material  number  in  the  li 
brary. 


0 


Code  to  specify  library  material. 


ID  AM 


1  will  compute  material  damage.  The  option  in  the  library  is 
IDAM  =  0,  which  will  not  compute  damage. 


I  FAIL 


DFRAC 


EFAIL 


=  0  will  not  allow  failure  of  the  material  when  the  damage  exceeds 
1.0,  but  rather  will  continue  to  accumulate  the  damage.  This  is 
the  option  which  is  in  the  library. 

1  will  allow  the  material  to  fracture  partially  when  the  damage 
exceeds  1.0.  Partial  fracture  causes  shear  an*-  tensile  failure,  so 
only  compressive  hydrostatic  pressure  capability  remains. 

-  Factor  by  which  library  fracture  strain  constants  <Dl,  D2,  EF 
MIN  —  defined  later)  are  multiplied.  If  left  blank  (DFRAC  =0) 
factor  will  be  set  to  DFRAC  =  1.0 

=  Equivalent  plastic  strain  (true)  which,  if  exceeded,  will  totally 
fail  the  element  such  that  it  produces  no  stresses  or  pressures.  If 
left  blank,  the  library  value  of  EFAIL  =  999.  will  govern. 


5  Material  Cards  for  Solids  Input  Data  j4I5,  5X,  F5.0,  A48/4<8F10.0)|  —  These 
five  cards  specify  all  the  material  constants  for  a  solid  material.  These  cards  will  su¬ 
persede  any  material  library  data  with  the  same  material  number,  MATL.  Only 
previously  undefined  variables  will  be  defined. 


DENSITY 

SPH  HEAT 

CONDUCT 

ALPHA 

TEMPI 

TROOM 

TMELT 

XI 

SHEAR  MOD 


Material  density  (mass/volume) 

Specific  heat  (work/mass/degree) 

Thermal  conductivity  (power/distance/degree) 

Volumetric  coefficient  of  thermal  expansion  (degree  ') 

Initial  temperature  of  the  material  (degree) 

Room  temperature  (degree) 

Melting  temperature  of  the  material  (degree) 

Extra  material  constant  stored  in  the  C9  array  (not  currently 
used) 

Shear  modulus  of  elasticity  (force  area) 


C1,C2,N,C3,  -  Constants  to  describe  the  material  strength,  o.  The  primary 

M.  C4,  SMAX  form  of  the  strength  equation  is 

a  =  |C1  s  C2 •  |1  +  C3 •  Cm *1  |1  -  T*mi  t  C4  P  ill 

Where  «  is  the  equivalent  plastic  strain,  t  =  «/<„  is  the 
dimensionless  strain  rate  for  tQ  =  1.0  s  ’,  T*  is  the  homologous 
temperature,  and  P  is  hydrostatic  pressure  (compression  is 
positive).  N  must  be  a  positive  number,  and  the  thermal  soften¬ 
ing  fraction,  Kq-  -  jl  —  T*M|,  is  set  to  Kq-  =1.0  when  M  =  0.  If 
SMAX  is  input  as  a  positive  number,  then  the  maximum 
strength  for  a  is  limited  to  SMAX.  If  left  blank  (SMAX  =  0.)  the 
strength  (a)  is  not  limited.  A  more  detailed  description  of  the 
strength  model  and  data  is  given  in  Reference  10. 

A  constant  flow  stress  can  be  obtained  by  setting  Cl  to  the  flow 
stress,  N  =  1.0,  and  C2  =  C3  =  C4  =  SMAX  =  M  =  0.  Cl,  C2 
and  SMAX  have  units  of  stress  (force/area)  and  the  others  are 
dimensionless. 

K1.K2,  K3  =  Cubic  coefficients  for  the  Mie-Gruneisen  equation  of  state 
(force/area) 

P  =  (Klg  +  K2p2  +  K3g3)  (1  —  fg/2)  +  rEs(l  +  g)  (2) 

where  g  =  p/p0  —  1  and  Es  is  internal  energy  per  initial  volume. 

F  =  Gruneisen  coefficient  for  Mie-Gruneisen  equation  of  state.  For 

the  average  pressure  options  (AVEP  =  2,4  in  the  Miscellaneous 
Card)  energy  effects  are  not  included,  so  it  is  necessary  to  set  I' 
=  0. 

CL  =  Linear  artificial  viscosity  coefficient  (CL  *  0.2). 

CQ  =  Quadratic  artificial  viscosity  coefficient  (CQ  =  4.0). 

PMIN  Maximum  hydrostatic  tension  allowed  iforce/area). 

CH  Hourglass  artificial  viscosity  coefficient  for  quad  elements  (CH 

*  0.02). 

1)1  D5  Constants  for  the  fracture  model  (Reference  llr 

iDl  *  D2expf,;i"'1  1  r  D4  •  Cm*  1  *  D5  •  T*  <3> 
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Where  f*  is  the  equivalent  strain  to  fracture  under  constant  con¬ 
ditions  of  the  dimensionless  strain  rate,  t  ,  homologous  tempera 
ture,  T  ,  and  the  pressure-stress  ratio,  a  ---  om/7T.  The  mean  nor¬ 
mal  stress  is  om  and  o  is  the  von  Mises  equivalent  stress.  Ex¬ 
pression  is  valid  for  a*  <  1.5.  Damage  is  computed  from  D  - 
and  fracture  is  allowed  to  occur  when  D  =  1.0. 

SPALL  Tensile  spall  stress  (negative  pressure)  at  which  fracture  can  oc¬ 

cur  (force/area). 

EFM1N  -  Minimum  fracture  strain  allowed.  For  o'  >  1.5,  varies  line¬ 

arly  from  (t  at  o*  =  1.5  to  EFMIN  at  am  -  SPALL. 

X2  =  Extra  material  constant  stored  in  the  CIO  array. 

Material  Cards  for  Explosives  from  Library  (215)  —  Similar  to  the  cards  for  the 

solid  materials  in  the  library  except  that  no  options  are  provided  for  fracture. 

3  Material  Cards  for  Explosives  Input  Data  (215,  20X,  A48/6F10.0/7F10.0)  — 

Only  new  variables  will  be  defined.  See  solid  material  definitions  for  other  vari¬ 
ables. 

ENERGY  Initial  internal  energy  in  explosive,  E(1  (energy/volume), 

DET  YEL  -  Detonation  velocity.  D  (distance/time). 

C1...C5  =  Constants  for  the  JWL  equation  of  state.  If  left  blank,  a  gamma 

law  equation  of  state  is  used  where 


For  gamma  law  the  pressure  is  determined  from 


P  -  (-> — 1)  E/V  (5) 

For  JWL  the  pressure  is  determined  from 

P  CI  -  (1-C5/C2  V)  •  exp(— C2  •  VI 
+■  CJ  ■  ( 1—-C5/C4  V)  ■  exp(— C4  •  V) 

*  C5  ■  EW  (6) 

where  E  is  internal  energy  per  initial  volume  and  V  V/V„  is 
the  relative  volume.  Cl  and  CJ  have  the  units  of  pressure 
(forcearea)  and  C2,  C4,  C5  are  dimensionless. 


XL  X2 


Extra  material  variables  stored  in  arrays  C6  and  C7 


Material  Cards  for  Crushable  Solids  from  Library  (215,  15X,  F5.0)  -  Similar  to 
the  cards  for  other  library  materials.  Total  failure  is  allowed  through  EFAIL,  but 
fracture  due  to  damage  is  not  allowed. 

4  Material  Cards  for  Crushable  Solids  Input  Data  |215,  15X,  F5.0, 
A48/3F10.0/2(8F10.0)|  —  These  four  cards  specify  the  material  constants  for  a 
crushable  solid  material  (Reference  12).  Only  new  variables  will  be  defined.  See 
previous  material  definitions  for  other  variables. 

Cl,  C4  =  Constants  to  describe  the  material  strength,  a 

a  -  Cl  +  C4  •  P  (7) 

SMAX  =  Maximum  strength  allowed  (force/area).  If  left  blank  (SMAX  = 

0.),  strength  ( a )  is  not  limited. 

PCRUSH,  =  Constants  to  describe  the  pressure,  P.  The  model,  and  specific 
UCRUSH,  data  for  concrete,  are  shown  in  Figure  15.  The  basic  model  can 

Kl,  K2,  K3,  also  be  used  for  other  crushable  solid  materials.  PCRUSH, 

KLOCK,  Kl,  K2,  K3  and  KLOCK  have  units  of  pressure  (force/areal. 

ULOCK  UCRUSH  and  ULOCK  are  dimensionless. 

XI  =  Extra  material  constant  stored  in  the  Dl  array.  End  material  in 

put  data  with  a  blank  card. 

Node  Geometry 

Node  geometry  data  are  required  for  the  projectile  nodes  and  the  target  nodes  These 
data  can  be  input  as  lines  of  nodes,  various  rod  shapes,  nose  shapes,  fiat  plates 
and/or  spheres.  The  input  data  are  summarized  in  Figure  3.  The  axial  coordinate,  z, 
is  positive  upward  and  the  radial  coordinate,  r,  is  positive  to  the  right.  One 
dimensional  geometry  (GEOM  =  3)  is  taken  along  the  z  axis  at  r  =  0. 

Line  of  Nodes  Description  Card  (415,  2X,  311,  5X,  5F10.0)  —  One  card  is  required 
for  each  line  of  nodes  to  be  generated.  No  additional  cards  are  needed.  The  nodes 
may  be  numbered  consecutively  or  incremented  by  INC  and  the  nodes  may  be  uni 
formly  or  variably  spaced.  Refer  to  Figure  8  for  more  details. 

1  Identification  number  for  line  of  nodes  geometry. 

N1  =  Number  of  the  first  node  of  the  line. 

NN  ~  Number  of  the  last  node  of  the  line.  Leave  blank  if  a  single  node 

is  to  be  generated. 


s  A.  J. 
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Node  number  increment  between  corresponding  nodes,  Leave 
blank  for  consecutive  numbering  or  if  a  single  node  is  to  be  gen¬ 
erated. 

=  Radial  restraint.  If  IR  =  1,  all  nodes  will  be  restrained  in  the  ra¬ 
dial  direction.  Leave  blank  for  no  restraint. 

=  Axial  restraint.  If  IZ  =  1,  all  nodes  will  be  restrained  in  the  ax¬ 
ial  direction.  Leave  blank  for  no  restraint. 

--  Theta  restraint.  If  IT  =  1,  all  nodes  will  be  restrained  in  the  0 
direction.  Leave  blank  for  no  restraint. 

=  R  coordinate  of  the  first  node. 

=  Z  coordinate  of  the  first  node. 


EXPAND 


R  coordinate  of  the  last  node.  Leave  blank  if  a  single  node  is  to 
be  generated. 

~  Z  coordinate  of  the  last  node.  Leave  blank  if  a  single  node  is  to 
be  generated. 

-  Factor  by  which  the  distance  between  nodes  is  multiplied  going 
from  the  first  node  to  the  last  node.  Leave  blank  for  uniform 
spacing. 


Rod  Node  Description  Card  (615,  3F10.0)  —  Two  cards  are  required  for  each  rod 
shape  to  be  generated.  Additional  cards  are  needed  for  certain  options.  The  first 
node  is  at  RMIN,  ZMAX,  and  the  nodes  are  numbered  across  each  layer  working 
down.  Radial  restraints  are  provided  when  the  inner  radius  lies  on  the  axis  of  sym 
metry  (r  —  0).  The  nodes  on  the  top  and  bottom  surfaces  of  the  rod  may  be  generated 
uniformly,  read  in  individually,  or  computed  by  analytic  functions.  Either  the  pri¬ 
mary  only  or  both  the  primary  and  secondary  (crossed  triangle)  nodes  may  be  gener¬ 
ated  in  the  rod  geometry.  (Refer  to  Figure  10  for  more  details). 

2  Identification  number  for  rod  nodes  geometry. 


NRING 


Number  of  the  first  node  of  the  rod  shape. 

Number  of  radial  rings  of  nodes  |not  including  the  secondary 
(crossed  triangle)  nodes  if  ICROS  ~  1|.  The  maximum  allowed  is 
50. 
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NPLN 


Number  of  horizontal  planes  of  nodes  jnot  including  the  secon¬ 
dary  (crossed  triangle)  nodes  if  ICROS  =  1|. 


r 


gives  uniform  radial  spacing  and  constant  z  coordinates  at  the 
top  and  bottom  of  the  rod. 


requires  all  r  and  z  coordinates  at  top  and  bottom  of  rod  to  be  in 
put  individually. 


IRAI)  -- 


< 


requires  r  coordinates  at  the  top  and  bottom  of  the  rod  to  be 
input  individually  and  z  coordinates  to  be  generated  with  an 
analytic  function. 


gives  uniform  radial  spacing  at  the  top  and  bottom  of  the  rod 
and  requires  z  coordinates  to  be  input  individually. 


gives  uniform  radial  spacing  at  the  top  and  bottom  of  the  rod 
and  requires  z  coordinates  to  be  generated  with  an  analytic  func¬ 
tion. 


ICROS 


Crossed  triangle  option.  If  ICROS  =  1,  both  the  primary  and  sec 
ondary  (crossed  triangle)  nodes  are  generated.  Leave  blank  to 
generate  the  primary  nodes  only. 


ZTOP 


Z  coordinate  of  the  top  surface  for  IRAD  =  0.  Leave  blank  for 
IRAD  =  1,2,  3,  4. 


ZBOT 


Z  coordinate  of  the  bottom  surface  for  IRAD  =  0.  Leave  blank  for 
IRAD  =  1,  2,  3,  4. 


EXPAND 


Factor  by  which  the  distance  between  corresponding  nodes  in 
the  vertical  direction  is  multiplied  going  from  top  to  bottom. 
Leave  blank  for  uniform  spacing  in  the  vertical  direction. 


Rod  Node 

ROTOP 

RITOP 

ROBOT 

RIBOT 


Radii  Card  for  IRAD  =  0  or  3  or  4  (4F10.0)  — 

=  Outer  radius  of  the  rod  top. 

-  Inner  radius  of  the  rod  top. 

Outer  radius  of  the  rod  bottom. 

Inner  radius  of  the  rod  bottom. 
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Rod  Node  Top  Radii  Card(s)  for  IRAD  =  1  or  2  (8F10.0)  — 

KT(  It ...  -  Radius  of  each  ring  of  nodes  at  the  top  of  the  rod.  One  or  more 

RT(NRING)  cards  as  required. 

Rod  Node  Bottom  Radii  Card(s)  for  IRAD  =  1  or  2  (8F10.0)  — 

RB(  1) ...  -  Radius  of  each  ring  of  nodes  at  the  bottom  of  the  rod.  One  or 

RBiNRING)  more  cards  as  required. 

Rod  Node  Top  Surface  Card(s)  for  IRAD  =  1  or  3  (8F10.0)  — 

ZT(  1 ) ...  =  Top  z  coordinate  of  each  ring  of  nodes.  One  or  more  cards  as 

ZT(NRING)  required. 

Rod  Node  Bottom  Surface  Card(s)for  IRAD  =  1  or  3  (8F10.0)  — 

ZB(1) ...  =  Bottom  z  coordinate  of  each  ring  of  nodes.  One  or  more  cards  as 

ZB(NRING)  required. 

Rod  Node  Top  Surface  Card  for  IRAD  =  2  or  4  (8F10.0)  — 

A,).  A!  ...  Ay  -  Coefficients  of  the  analytical  function  describing  the  top  surface 

|zlop  =  A0  +  Ajr  +  ...  +  A6r6  +  A7  (1  -  cos  0)| 

Rod  Node  Bottom  Surface  Card  for  IRAD  =  2  or  4  (8F10.0)  — 

B0.Bj  ...  B7  =  Coefficients  of  the  analytical  function  describing  the  bottom  sur¬ 
face  |z|)((t  -  B0  +  B]r  +  ...  +  B6r6  +  B7  0  -  cos  0S| 

Note:  If  it  is  not  possible  to  describe  the  node  geometry  of  the  rod  with  a  single 
shape,  it  is  possible  to  use  multiple  shapes  to  form  a  single  rod.  The  nodes  must  he 
numbered  consecutively,  and  the  radii  and  the  number  of  rings  of  nodes  must  be  the 
same  for  the  individual  rod  shapes  at  their  interface.  Also  ZTOP  and  ZBOT  for  the 
adjoining  rods  should  not  be  identical.  ZTOP  for  the  lower  rod  should  be  less  than 
ZBOT  for  the  upper  rod  by  the  desired  nodal  spacing  in  the  z  direction.  If  the  crossed 
triangle  option  is  being  used,  the  secondary  nodes  between  the  two  shapes  must  be 
generated  individually  or  with  a  line  of  nodes. 

Nose  Node  Description  Card  (615,  3F10.0)  —  One  card  is  required  for  each  nose 
shape  to  be  generated.  Additional  cards  are  needed  for  certain  options.  The  first 
node  <NIi  is  the  r  0,  z  -  ZMIN  node  generated  by  the  rod  generator.  The  nose  gen¬ 
erator  numbers  the  nodes  consecutively  in  the  clockwise  direction,  starting  at  the 
central  node  and  working  outward,  ring  by  ring.  The  nose  shapes  are  always  getter 
ated  pointing  downwards.  The  number  of  rings  of  nodes  for  the  nose  must  be  one  less 


than  the  number  of  rings  of  nodes  for  the  rod.  Either  primary  only  or  primary  and 

secondary  (crossed  triangle)  nodes  may  he  generated  in  the  nose  geometry  Refer  to 

Figure  1 1  for  more  details. 

3  —  Identification  number  for  nose  nodes  geometry. 

N1  -  Number  of  the  first  node  at  the  rod-nose  interface. 

NRING  =  Number  of  rings  |not  including  the  secondary  (crossed  triangle) 

nodes  if  ICROS  =  1|.  The  maximum  allowed  is  25. 

INOSE  -  Nose  option.  If  INOSE  =  1,  a  rounded  nose  is  generated.  If  the 

length  of  the  nose  is  equal  to  the  radius,  a  hemispherical  nose  is 
generated.  If  INOSE  =  2,  a  conical  nose  is  generated.  If  INOSE 
=  3,  a  tangent  ogival  nose  is  generated.  The  length  of  the  ogival 
nose  must  be  greater  than  or  equal  to  the  radius  at  the  nose-rod 
interface. 

IRAD  =  Radius  option.  If  IRAD  =  1,  individual  radii  and  z  coordinates 

are  input  for  each  ring  of  nodes.  Leave  blank  (IRAD  =  0)  for  uni 
form  spacing. 

ICROS  -  Crossed  triangle  option.  If  ICROS  1,  both  the  primary  and  see 

ondary  (crossed  triangle)  nodes  are  generated.  Leave  blank  to 
generate  the  primary  nodes  only. 

RTOP  -  R  coordinate  of  the  outer  ring  of  nodes  at  the  rod-nose  interface. 

Leave  blank  for  IRAD  =  1. 

ZTOP  -  Z  coordinate  at  the  rod-nose  interface  and  thus  equal  to  ZBOT 

from  the  rod  geometry. 

ZBOT  -  Z  coordinate  at  the  tip  of  the  nose.  Leave  blank  for  IRAD  =  1 . 

Nose  Node  Top  Radii  Card(s)  for  IRAD  =  1  (8F10.0)  — 

RTili...  ~  R  coordinate  of  each  ring  of  nodes  at  the  rod-nose  interface  (not 

RT  (NRING)  including  the  central  node  which  is  on  the  axis  of  symmetry). 

One  or  more  cards  as  required. 

Nose  Node  ZMIN  Card(s)  for  IRAD  =  1  (8F10.0)  - 

ZBi  1 1 ...  Z  coordinate  at  centerline  of  each  ring  of  nodes  at  the  axis  of 

ZB  (NRING)  symmetry  (not  including  the  central  node  which  is  at  ZTOP). 

One  or  more  cards  as  required. 


Flat  Plate  Description  Card  (615,  4F10.0)  —  Two  cards  are  required  for  each  flat 
plate  shape  to  be  generated.  No  additional  cards  are  needed.  The  first  node  is  at  the 
RMIN,  ZMAX  corner  of  the  plate  and  the  nodes  are  numbered  across  the  plate  work¬ 
ing  down.  Flat  plates  have  horizontal  tops  and  bottoms  and  vertical  sides.  They 
made  be  joined  together  side  to  side,  with  or  without  crossed  triangles,  if  the  join  op¬ 
tion  is  used  and  the  node  number  increment  is  equal  to  the  total  number  of  nodes 
(primary  and  secondary)  in  the  radial  direction.  Regions  with  variable  nodal  spacing 
may  he  included  at  the  RMAX  end  and/or  the  ZMIN  end.  Refer  to  Figure  12  for  more 
details. 

4  =  Identification  number  for  flat  plate  nodes  geometry. 

N 1  -  Number  of  the  first  node  of  the  plate. 

INC  =  Node  number  increment  between  corresponding  nodes  in  the 

vertical  direction.  Leave  blank  if  the  join  option  is  not  being 
used. 

IRFIX  Radial  restraint  option.  If  IRFIX  =  1,  all  nodes  on  the  axis  of 

symmetry  are  restrained  radially.  Leave  blank  for  no  restraint. 

IJOIN  =  Join  option  or  crossed  triangle  geometry  (ICROS  =  1).  If  IJOIN 

=  I,  the  generation  of  the  first  column  of  primary  nodes  is  sup¬ 
pressed  so  that  plates  may  be  joined  together  side  to  side.  Leave 
blank  if  only  primary  nodes  are  generated. 

ICROS  -  Crossed  triangle  option.  IF  ICROS  =  1,  both  the  primary  and 

secondary  (crossed  triangle)  nodes  are  generated.  Leave  blank  to 
generate  the  primary  nodes  only. 

RMAX  =  Maximum  r  coordinate  of  the  plate. 

RMIN  =  Minimum  r  coordinate  of  the  plate. 

ZMAX  Maximum  z  coordinate  of  the  plate. 

ZMIN  Minimum  z  coordinate  of  the  plate. 

Hal  Plate  Expansion  Card  (415,  10X,  4F10.0)  — 

NR  Number  of  nodes  in  the  r  direction  not  including  the  secondary 

(crossed  triangle)  nodes  if  ICROS  lj.  The  maximum  allowed  is 
50. 

NZ  Number  of  nodes  in  the  z  direction  inot  including  the  secondary 

(cross  triangle)  nodes  if  ICROS  1  , 


NA 


NREND 

NZEND 

RPART 

R  EXPAND 

ZPART 

Z-EXPAND 


-  Number  of  nodes  in  the  r  variable  node  spacing  section.  The 
node  at  the  division  between  the  uniform  and  the  variable  spac 
ing  sections  is  included  in  this  number.  Leave  blank  for  uniform 
spacing  in  the  r  direction  and  set  equal  to  NR  for  variable  spac¬ 
ing  only. 

=  Number  of  nodes  in  the  z  variable  node  spacing  section.  The 
node  at  the  division  between  the  uniform  and  the  variable  spac¬ 
ing  sections  is  included  in  this  number.  Leave  blank  for  uniform 
spacing  in  the  z  direction  and  set  equal  to  NZ  for  variable  spac¬ 
ing  only. 

=  Fractional  part  of  the  radial  length  occupied  by  the  variable 
spacing.  Leave  blank  for  uniform  spacing  in  the  radial  direction. 

=  Factor  by  which  the  radial  distance  between  nodes  is  multiplied 
in  the  variable  spacing  section  moving  outward.  Leave  blank  for 
uniform  spacing  in  the  r  direction. 

=  Fractional  part  of  the  axial  length  occupied  by  the  variable  spac 
ing.  Leave  blank  for  uniform  spacing  in  the  axial  direction. 

=  Factor  by  which  the  z  distance  between  nodes  is  multiplied  in 
the  z  variable  spacing  section  moving  downward.  Leave  blank 
for  uniform  spacing  in  the  z  direction. 


Sphere  Node  Description  Card  (415,  10X,  2F10.0)  —  One  card  is  required  for  each 
sphere  shape  to  be  generated.  No  additional  cards  are  needed.  The  first  node  is  the 
central  node  and  the  nodes  are  numbered  consecutively  in  the  clockwise  direction, 
starting  at  the  central  node  and  working  outwards  ring  by  ring.  Either  primary  only 
or  both  primary  and  secondary  (crossed  triangle)  nodes  may  be  generated  in  the 
sphere  geometry.  Refer  to  Figure  13  for  more  details. 


5 

Nl 

NRING 

ICROS 

ZTOP 

ZROT 


Identification  number  for  sphere  nodes  geometry. 

Number  of  the  first  node  of  the  sphere. 

Number  of  rings  of  nodes  |not  including  the  secondary  (crossed 
triangle)  node  if  ICROS  =  1|.  The  maximum  allowed  is  25. 

Crossed  triangle  option.  If  ICROS  =  1,  both  the  primary  and  .sec¬ 
ondary  (crossed  triangle)  nodes  are  generated.  Leave  blank  to 
generate  the  primary  nodes  only. 

Top  z  coordinate  of  the  sphere 

Bottom  z  coordinate  of  the  sphere 
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NABOR  Node  Description  Card  (315,  15X,  5F10.0)  —  One  card  is  required  for 
each  NABOR  node  shape  to  be  generated.  A  description  of  the  NABOR  node  for 
mulation  will  appear  in  the  final  report  for  this  contract.  The  NABOR  node  option  is 
available  for  axisymmetric  or  plane  strain  geometry  (GEOM  =  1  or  2).  However, 
there  is  no  allowance  for  rotational  displacements  in  the  6  direction  (about  the  z  axis 
of  symmetry).  Only  solid  materials  or  crushable  solid  materials  (no  explosive  materi¬ 
als)  can  be  used  as  NABOR  nodes,  and  the  two  different  material  types  cannot  be  in 
contact.  NABOR  nodes  may  be  used  exclusively,  or  they  may  be  used  in  conjunction 
with  the  standard  nodes  and  elements.  The  mass  of  NABOR  nodes  is  dependent  only 
on  the  density  of  the  material,  the  radial  coordinate  (for  GEOM  =  1)  and  the 
diameter  of  the  node.  They  cannot  accept  additional  mass  from  attached  elements  or 
input  concentrated  masses.  The  initial  diameters  of  all  NABOR  nodes  in  a  specific 
problem  must  be  identical. 

Figure  14  shows  initial  geometry  for  NABOR  nodes.  They  are  arranged  in  a  hex¬ 
agonal  manner.  The  first  node,  Nl,  is  upper-left,  adjusted  as  shown.  Nodes  are 
numbered  consecutively,  from  left  to  right  in  horizontal  rows,  moving  downward. 
The  even-numbered  rows  have  one  less  node  than  the  odd-numbered  rows.  Other  ori¬ 
entations  may  be  achieved  by  using  the  Rotate  option  in  the  Penetrator  and/or 
Target  Scale/Shift/Rotate  Cards. 

6  =  Identification  number  for  NABOR  nodes  geometry. 

Nl  =  Number  of  first  node. 

MATL  =  Material  number  for  the  NABOR  nodes. 

RMAX  =  Maximum  r  coordinate  of  the  NABOR  node  shape. 

RMIN  Minimum  r  coordinate  of  the  NABOR  node  shape. 

ZMAX  -  Maximum  z  coordinate  of  the  NABOR  node  shape. 

ZMIN  =  Minimum  z  coordinate  of  the  NABOR  node  shape. 

NODE  DIA  =  NABOR  node  diameter.  Must  be  identical  for  all  NABOR  nodes. 

Can  be  left  blank  for  subsequent  NABOR  node  shapes. 

Element  Geometry 

The  element  data  are  required  to  be  consistent  with  the  node  data  for  the  projectile 
and  the  target.  Thus,  a  series  of  composite  elements,  rod  elements,  nose  elements, 
Hal  plate  elements,  and  sphere  elements  may  be  created.  The  element  data  for  these 
shapes  are  entered  individually  in  the  locations  identified  in  Figure  1.  There  is  no 
limit  to  the  number  of  special  shapes  that  may  be  used  in  the  projectile  or  the  target. 


V*  \’a  .  , 


If  any  quad  elements  are  used,  there  is  no  allowance  for  rotational  displacements  in 
the  6  direction  (about  the  z  axis  of  symmetry).  The  input  data  are  summarized  in 
Figure  4. 


Series  of  Composite  Elements  Description  Card  (1015)  —  One  card  is  required 
for  each  series  of  composite  elements  to  be  generated.  No  additional  cards  art- 
needed.  Standard  or  crossed  triangle  elements  are  generated  depending  on  the 
number  of  nodes  input.  One-dimensional  elements  and  two-dimensional  quad  ele¬ 
ments  can  also  be  generated.  The  elements  are  numbered  consecutively  and  only  one 
material  may  be  used  per  series.  For  triangle  and  quad  elements,  the  nodes  are  in 
put  in  a  counterclockwise  manner.  For  one-dimensional  geometry,  node  N1  must 
have  a  higher  z  coordinate  than  node  N2.  Refer  to  Figure  9  for  more  details. 


NCOMP 
MATL 
N 1 


=  Identification  number  for  series  of  elements  geometry. 
=  Number  of  composite  elements  to  be  generated. 

=  Material  number  for  the  series  of  elements. 

=  Number  of  the  first  node. 

=  Number  of  the  second  node. 


=  Number  of  the  third  node.  Leave  blank  for  one-dimensional 
geometry  (GEOM  =  3). 


=  Number  of  the  fourth  node.  Leave  blank  if  a  single  triangular 
element  is  to  be  generated  per  composite  element. 


Number  of  the  fifth  node  for  crossed  triangles.  Leave  blank  if 
one  or  two  triangular  elements  are  to  be  generated  per  com¬ 
posite  element.  Set  N5  =  —  1  for  quad  elements. 


=  Node  number  increment  between  corresponding  nodes  in  each 
composite  element.  Leave  blank  if  the  nodes  are  numbered  con¬ 
secutively. 


Number  of  the  first  element  in  the  series  of  elements 


Example:  To  generate  the  elements  shown  at  the  bottom  of  Figure  9,  the  following 
data  must  be  used. 


N1 

N2 

N3 

N4 

106 

101 

102 

107 
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Rod  Element  Description  Card  (715)  —  One  card  is  required  for  each  rod  shape  to 
lie  generated.  Additional  cards  are  needed  for  the  material  option.  The  elements  are 
numbered  consecutively  across  the  rod,  working  down  layer  by  layer.  Standard,  al¬ 
ternating  diagonal,  crossed  triangle  elements,  or  quad  elements,  may  he  generated. 
Refer  to  Figure  10  for  more  details. 

2  Identification  number  for  rod  elements  geometry. 


MAIL 


Number  of  the  first  node  of  the  rod  shape. 

Material  number  for  the  rod  elements.  Leave  blank  if  each  ver 
tieal  column  of  composite  elements  is  to  have  a  different  materi¬ 
al  number. 


NRCOL 


NZLAY 


IDIAG 


-  Number  of  vertical  columns  of  composite  elements.  Note  that 
there  is  one  less  column  of  elements  when  compared  to  the 
number  of  rings  of  nodes.  The  maximum  allowed  is  50. 

Number  of  horizontal  rows  of  c  omposite  elements. 

=  Diagonal  option.  If  IDIAG  5,  it  must  be  used  with  the  crossed 
triangle  option  for  nodal  geometry  (ICROS  -  D.  Refer  to  Figure 
10. 


LI  Number  of  the  first  element  in  the  rod  shape. 

Material  Card(s)  for  MATL  =  0  (1615)  — 

Mill...  Material  number  for  each  vertical  column  of  composite  ele- 

MiNRCOL)  ments.  Numbered  from  inner  column  to  outer  column.  One  or 

more  cards  as  required. 

Nose  Element  Description  Card  (415,  5X,  215)  —  One  card  is  required  for  each 
nose  shape  to  be  generated.  Additional  cards  are  needed  for  the  material  upturn.  The 
elements  are  numbered  consecutively  in  the  clockwise  direction,  starting  at  the  first 
node  and  working  outwards  ring  by  ring.  Standard  or  crossed  triangle  eh  u  tits,  or 
quad  elements,  may  be  generated.  Refer  to  Figure  1 1  for  more  details. 

•'!  Identification  number  for  nose  elements  geometry 

Nl  Number  of  the  first  node-  of  the  nose  shape  at  nose  rod  interface. 


MAIL 


Material  number  for  tile  nose  elements.  Leave  blank  if  each  ra 
dial  ring  of  composite  elements  is  to  have  a  different  material 
number. 


.N  N  --  \ 


NRCOL  =  Number  of  radial  rings  of  composite  elements.  The  maximum  al¬ 

lowed  is  25. 

IDIAG  -  Diagonal  option.  If  IDIAG  =  5,  it  must  be  used  with  the  crossed 

triangle  geometry  for  nodal  geometry  (ICROS  =1).  Set  IDIAG  = 
0  for  standard  triangular  elements  and  IDIAG  —  6  for  quad  ele¬ 
ments. 

LI  =  Number  of  the  first  element  in  the  nose  shape. 

Material  Card(s)  for  MATL  =  0  (1615)  — 

M(l)...  =  Material  number  for  each  radial  ring  of  composite  elements.  One 

M(NRCOL)  or  more  cards  as  required. 

Flat  Plate  Element  Description  Card  (715)  —  One  card  is  required  for  each  flat 
plate  shape  to  be  generated.  No  additional  cards  are  needed.  The  elements  are 
numbered  consecutively  across  the  plate,  working  down  layer  by  layer.  Standard,  al¬ 
ternating  diagonal,  or  crossed  triangle  elements  may  be  generated.  Refer  to  Figure 
12  for  more  details. 

4  =  Identification  number  for  flat  plate  element  geometry. 

N1  =  Number  of  the  first  node  of  the  flat  plate  shape. 

MATL  =  Material  number  for  the  flat  plate  elements. 

NRCOL  =  Number  of  vertical  columns  of  composite  elements. 

NZLAY  -  Number  of  horizontal  layers  of  composite  elements. 

IDIAG  =  Diagonal  option.  If  IDIAG  =  5,  it  must  be  used  with  the  crossed 

triangle  option  for  nodal  geometry  (ICROS  =1).  Refer  to  Figure 

12. 

LI  =  Number  of  the  first  element  of  the  flat  plate  shape. 

Sphere  Element  Description  Card  (415, 5X, 215)  —  One  card  is  required  for  each 
sphere  shape  to  be  generated.  No  additional  cards  are  needed.  The  elements  are 
numbered  consecutively  in  the  clockwise  direction,  starting  at  the  central  node  and 
working  outwards  ring  by  ring.  Standard  or  crossed  triangle  elements,  or  quad  ele 
ments  may  be  generated.  Refer  to  Figure  13  for  more  details. 

5  =  Identification  number  for  sphere  elements  geometry. 

N1  =  Number  of  the  first  node  of  the  sphere  shape. 


MATL 


Materia)  number  for  the  sphere  elements. 


NRCOL  =  Number  of  radial  ring's  of  composite  elements.  The  maximum  al¬ 

lowed  is  25. 

1 1)1  AC  1  =  Diagonal  option.  IF  1DIAG  5.  it  must  be  used  with  the  crossed 

triangle  option  for  nodal  geometry  (1CROS  -  It  Set  IDIAG  0 
for  standard  triangular  elements  and  IDIAG  -  6  for  quad  ele¬ 
ments. 

Ll  =  Number  of  the  first  element  of  the  sphere  shape. 

Material  Card(s)  for  MATL  =  0  (1615)  — 

-  Material  number  for  each  radial  ring  of  composite  elements.  One 
Mi  NRCOL i  or  more  cards  as  required. 

INPUT  DATA  FOR  THE  MAIN  ROUTINE 

The  function  of  the  Main  Routine  is  to  perform  the  computations.  It  may  be  used  in 
conjunction  with  the  Preprocessor,  or  it  can  read  initial  conditions  from  the  restart 
tape  which  has  been  previously  generated  from  a  Preprocessor  run  or  another  Main 
Routine  run.  The  descriptions  which  follow  are  for  the  data  in  Figure  5.  Consistent 
units  must  be  used  and  the  unit  of  the  time  must  be  in  seconds. 

Restart  Description  Card  (215,  A 60)  —  This  card  is  used  only  for  restart  runs.  If 
the  Main  Routine  is  run  in  conjunction  with  the  Preprocessor  (TYPE  -  2  on 
Preprocessor  Miscellaneous  Card),  then  this  card  is  omitted. 

CASE  =  Case  number  for  run. 

.)  =-  Code  to  indicate  restart  run. 

DESCRIPTION  ---  Description  of  problem  provided  by  the  user. 

Time  Integration  Card  (2I5.7F10.0)  — 

r  Cycle  number  at  which  the  run  begins.  The  cycle  numbers  for 
which  restart  files  are  written  are  given  in  the  printed  output  of 
the  previous  run  (Preprocessor  or  Main  Routine).  CYCLE  =  0  if 
the  Main  Routine  is  used  in  conjunction  with  the  Preprocessor 


CYCLE 


0  will  not  allow  problem  size  to  be  changed. 

1  will  allow  the  problem  size  to  be  reduced  at  a  specified  time. 

(  HANG  <  2  will  allow  the  problem  size  to  be  reduced  at  one  time,  and  then 

increased  at  another  time.  The  expansions  and  reductions  can  be 
performed  in  either  order.  If  identical  times  are  designated  for 
the  expansion  and  reduction,  the  reduction  will  be  performed 
first. 

3  will  allow  the  problem  size  to  be  expanded  at  a  specific  time. 

DTMAX  =  Maximum  integration  time  increment  which  will  be  used  for  the 

equations  of  motion. 

DTMIN  =  Minimum  integration  time  increment  allowed.  If  exceeded,  the 

results  will  be  written  onto  the  restart  tape  and  the  run  will 
stop. 

SSF  =  Fraction  of  the  sound  speed  transmit  time  used  for  the  integra¬ 

tion  time  increment.  Must  be  less  than  1.0. 

'I'M AX  Maximum  time  the  problem  is  allowed  to  run.  This  time  refers 

to  the  dynamic  response  of  the  system,  not  the  central  processor 
time  (CPMAX)  described  next.  The  results  at  time  =  TMAX  are 
written  onto  the  restart  tape,  and  the  run  is  discontinued 

CPMAX  =  Central  processor  time  at  which  the  results  will  be  written  onto 

the  restart  tape  and  the  run  will  stop.  The  time  units  for  this  in 
put  can  be  seconds,  minutes  or  hours.  It  should  coincide  with  the 
units  the  specific  computer  uses  to  measure  central  processor 
time. 

KMAX  Upper  limit  for  total  kinetic  energy.  This  is  used  for  numerical 

instability  checks.  The  run  will  stop  if  the  kinetic  energy  ex 
ceeds  EMAX.  If  left  blank,  EMAX  will  automatically  be  set  to 
1.5  times  the  initial  total  energy. 

VNKKF  Reference  closing  velocity  for  NABOR  search  routines.  It  should 

be  0.2  to  1.0  times  the  impact  velocity  (or  relative  closing  voloc 
ity)  in  the  problem.  Lower  values  will  cause  less  searching  but 
may  allow  two  NABOR  nodes  to  become  too  close  before  being 
detected.  Higher  values  will  cause  more  searching  but  will  also 
ensure  that  nodes  do  not  become  too  close  before  being  detected. 


Change  Drop  Card  for  CHANG  =  1  or  2  (F10.0,  1315)  —  This  card  is  used  only  if 
changes  are  made  which  reduce  the  size  of  the  problem.  The  portions  of  the  problem 
which  remain  are  those  which  were  input  first.  Common  uses  are  to  drop  the  ex¬ 
plosive  gases  after  a  liner  has  been  accelerated,  or  to  drop  the  target  after  a 
projectile  has  perforated  the  target.  For  CHANG  =  1,  the  problem  size  is  reduced 
only  and  there  are  no  expansions.  For  CHANG  =  2,  there  will  be  problem  ex¬ 
pansions  as  well  as  problem  reductions. 


TDROP 


NODE 


PNODE 


PEEK 


NMAS 


NRIG 


NCHNK 


NSLID 


NPEOT 


EPEOT 


NFAIE 


PNAB 


=  Time  at  which  the  change  (problem  size  reduction)  occurs. 

-  Total  number  of  nodes  which  remain  in  the  revised  problem. 

Number  of  projectile  nodes  which  remain  in  the  revised  problem. 

Total  number  of  elements  which  remain  in  the  revised  problem 

■-  Number  of  projectile  elements  which  remain  in  the  revised  prob¬ 
lem. 

--  Number  of  concentrated  masses  which  remain  in  the  revised 
problem. 

-  Number  of  rigid  systems  of  nodes  which  remain  in  the  revised 
problem. 

Number  of  subsystems  of  chunks  of  elements  which  remain  in 
the  revised  problem. 

-  Number  of  sliding  surfaces  which  remain  in  the  revised  prob¬ 
lem. 

Number  of  nodes,  for  which  time-history  data  are  written,  which 
remain  in  the  revised  problem. 

Number  of  elements,  for  which  time  history  data  are  written, 
which  remain  in  the  revised  problem. 

Number  of  elements  which  will  be  designated  to  fail  totally 
This  type  of  failure  sets  all  stresses  in  the  element  to  zero.  It  es¬ 
sentially  makes  the  element  disappear  except  that  mass  is  re 
tained  at  the  nodes. 

=  Total  number  of  NABOR  nodes  which  remain  in  the  revised 
problem. 

Number  of  projectile  NABOR  nodes  which  remain  in  the  revised 
problem. 


-V.y.v.y/,- 
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Designated  Element  Failure  Card  (1615)  —  This  card  is  used  only  if  there  are  ele 
ments  to  be  totally  failed  (NFAIL  >0). 

EF1  ...  EFN  =  Elements  to  be  totally  failed  in  the  revised  problem 

Change  Add  Card  for  CHANG  =  2  or  3  (F10.0,  1115)  —  This  card  is  used  only  if 
there  are  additions  to  be  made  to  the  problem.  These  additions  can  be  made  in  con 
junction  with  previous  problem  reductions  (CHANG  =  2 1  or  they  can  be  made  to  the 
original  problem  (CHANG  =  3). 

TADD  =  Time  at  which  the  change  (problem  size  expansion)  occurs. 

NMAS  =  Number  of  concentrated  masses  in  the  expanded  problem.  If  this 

is  greater  than  the  existing  number  of  concentrated  masses,  the 
program  will  read  input  data  for  the  additional  masses.  As  an 
example  for  CHANG  =  2,  consider  the  case  of  NMAS  :  3  in  the 
initial  geometry  generated  by  the  preprocessor,  NMAS  =  2  in 
the  Change  Drop  Card  and  NMAS  =  6  here  in  the  Change  Add 
Card. 

The  initial  geometry  would  read  three  concentrated  masses. 
These  would  be  included  in  the  problem  until  time  =  TDROP,  at 
which  time  only  the  first  two  concentrated  masses  would  be  in 
eluded  (the  third  would  be  eliminated).  Now,  because  the  ex 
panded  problem  will  have  six  concentrated  masses,  the  program 
needs  to  read  data  for  an  additional  four  masses  (6  expanded  2 
existing  =  4  required).  These  masses  will  come  into  effect  at 
time  =  TADD. 

NRIG  =  Number  of  rigid  systems  of  nodes  in  the  expanded  problem  If 

this  is  greater  than  the  existing  number  of  rigid  systems  of 
nodes,  the  program  will  read  input  data  for  the  additional  ays 
terns  of  nodes. 

NCHNK  -  Number  of  subsystems  of  chunks  of  elements  in  the  expanded 

problem.  If  this  is  greater  than  the  existing  number  of 
subsystems  of  chunks  of  elements,  the  program  will  read  input 
data  for  the  additional  chunks  of  elements. 

NSLID  =  Number  of  sliding  surfaces  in  the  expanded  problem.  If  this  is 

greater  than  the  existing  numbers  of  sliding  surfaces,  the  pro¬ 
gram  will  read  input  data  for  the  additional  sliding  surfaces. 

IZR  =  1  gives  a  rigid  frictionless  surface  on  the  positive  side  of  the 

plane  described  by  z  =  0.  If  the  equations  of  motion  cause  a  node 
to  have  a  negative  z  coordinate,  the  z  coordinate  and  velocity  are 
set  to  zero.  If  IZR  -  0  this  option  is  not  used. 
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SPLIT 


0  or  1.  See  description  given  for  Miscellaneous  Card  in  the 
Preprocessor. 


NTOP,  NBOT,  =  Node  geometry  descriptions  if  internal  loads  are  desired.  See  de 
NRING  scription  given  for  Miscellaneous  Card  in  the  Preprocesor.  Musi 

be  input  here  even  if  input  in  initial  geometry. 


NRST 


IPRNT 


-  Number  of  groups  of  nodes  to  have  restraints  redefined. 

=  1  will  print  individual  node  and  element  geometry  of  the  ex 
panded  problem.  Will  not  print  data  for  IPRNT  =  0. 


Target  Scale/Shift/Rotate  Card  for  CHANG  =  2  or  3  (F10.0)  —  This  card  is  the 
same  as  used  in  Preprocessor.  Applies  to  all  target  nodes  in  the  expanded  problem. 
Must  be  included  even  if  there  are  no  additional  nodes  to  be  input. 


r\ 


CM 


1 


Additional  Target  Node  Data  for  CHANG  =  2  or  3  —  These  cards  are  required  to 
define  the  additional  nodes  (for  target  only)  for  the  expanded  problem.  Same  format 
as  used  in  Preprocessor.  End  with  a  blank  card.  Include  blank  card  even  if  there  are 
no  additional  target  nodes. 

Additional  Target  Element  Data  for  CHANG  =  2  or  3  —  These  cards  are  re 
quired  to  define  the  additional  elements  (for  target  only)  for  the  expanded  problem. 
Same  format  as  used  in  Preprocessor.  End  with  a  blank  card.  Include  blank  card 
even  if  there  are  no  additional  target  elements. 

NMASS.  NRIG,  NCHNK,  NSLID,  NRST  data  read  as  required.  Same  format  as  used 
in  Preprocessor. 

Target  Velocity  Card  for  CHANG  =  2  or  3  (30X,  3F10.0)  —  This  card  defines  the 
nodal  velocities  of  all  target  nodes  in  the  expanded  problem.  See  descriptions  for  Ve¬ 
locity/Detonation  Card  in  Preprocessor.  Must  be  included  even  if  there  are  no  target 
nodes  or  if  no  additional  target  nodes  were  input.  If  there  are  no  target  velocities,  do 
not  leave  the  card  blank.  Instead,  input  velocities  as  0.  Must  be  followed  by  a  blank 
card. 

Plot  Card  <4I5,4F10.0)  - 


PCYCL 


Same  restart  cycle  as  specified  on  the  Time  Integration  Card. 
0  will  not  write  the  system  data  on  the  plot  tape. 

1  will  write  all  the  system  data  on  the  plot  tape. 


• .v 
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N* 


NPLOT 


LPLOT 


DT  SYS 


Number  of  nodes  for  which  data  will  be  written  on  the  plot  tape. 
The  individual  nodes  are  specified  on  the  Designated  Nodes 
Cards.  Program  is  currently  dimensioned  for  a  maximum  of  20 
nodes. 

Number  of  elements  for  which  data  will  be  written  on  the  plot 
tape.  The  individual  elements  are  specified  on  the  Designated 
Elements  Cards.  Program  is  currently  dimensioned  for  a  max 
imum  of  20  elements. 

Time  increment  at  which  the  system  data  are  written  on  the  plot 
tape.  These  quantities  do  not  vary  as  rapidly  as  do  the  individ 
ual  node  and  element  data  so  a  larger  time  increment  can  be 
used. 


r~i 
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TSYS  =  Time  at  which  the  first  system  data  are  written  on  the  plot  tape. 

If  left  blank,  the  time  at  the  beginning  of  the  Main  Routine  run 
will  be  used. 

DT  NODE  =  Time  increment  at  which  the  individual  node  and  element  data 
are  written  on  the  plot  tape.  These  quantities  vary  more  rapidly 
than  the  system  data  so  a  smaller  time  increment  can  be  used. 

TNODE  =  Time  at  which  the  first  individual  node  and  element  data  are 

written  on  the  plot  tape. 

Designated  Nodes  Card  (1615)  —  This  card  is  used  only  if  there  are  node  data  to  be 

written  on  the  plot  tape  (NPLOT >0  on  the  plot  card). 

N1...NN  =  Individual  node  numbers  for  which  data  will  be  written  on  the 

plot  tape.  Must  be  in  ascending  order. 

Designated  Elements  Card  (1615)  —  This  card  is  used  only  if  there  are  element 

data  to  be  written  on  the  plot  tape  (LPLOT  >0  on  the  plot  card). 


E1...EN 


Individual  element  numbers  for  which  data  will  be  written  on 
the  plot  tape.  Must  be  in  ascending  order. 


Data  Output  Cards  (3F10.0,6I5)  —  These  cards  are  used  to  specify  various  forms  of 
output  data  at  selected  times,  and  the  last  card  must  be  for  a  time  greater  than 
TMAX  even  though  output  will  not  be  provided  for  that  specific  time.  Recall  that 
output  is  automatically  provided  at  TMAX,  and  a  data  output  card  need  not  he 
provided  for  this  time.  End  run  with  a  blank  card. 


TIME 


Time  at  which  output  data  will  be  provided. 
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EC HECK 


Code  which  governs  the  printed  output.  The  following  options 
are  provided: 


RDAMP 


SAVE 


LOAD 


CHUNK 


BURN 


NDATA 


PLOT 


v  v -j 


1.  If  ECHECK  is  greater  than  1000,  the  individual  node  data 
and  element  data  will  not  be  printed.  Only  system  data 
such  as  eg  positions,  momenta,  energies  and  average  veloc¬ 
ities  are  provided  for  the  projectile,  target  and  entire  sys¬ 
tem. 

2.  If  ECHECK  is  less  than  1000,  the  system  data  and  individ¬ 
ual  node  data  will  be  printed.  Individual  solid  element  data 
will  be  printed  for  all  elements  which  have  an  equivalent 
plastic  strain  equal  to  or  greater  than  ECHECK.  For  exam 
pie,  if  ECHECK  =  0.5,  all  elements  with  equivalent  plastic 
strains  equal  to  or  greater  than  0.5  will  have  data  printed. 
If  ECHECK  =  999,  no  element  data  will  be  printed. 

Radial  damping  constant,  Cp  in  equation  (38)  of  Reference  1 
This  damping  acts  until  the  time  specified  in  the  following  Data 
Output  Card.  RDAMP  =  0  for  GEOM=2,3,4  and  axisymmetric 
geometry  when  no  spin  is  present.  For  problems  involving  spin. 
Reference  1  should  be  consulted. 

1  will  write  results  on  restart  tape  for  possible  restart  runs  or 
state  plots.  Will  not  write  results  if  left  blank. 

1  will  compute  and  print  internal  loads  in  a  cylindrical 
projectile.  This  option  can  only  be  used  if  NTOP,  NBOT  and 
NRING  are  properly  defined  in  the  Miscellaneous  Card  of  the 
Preprocessor.  Will  not  compute  if  left  blank. 

1  will  compute  and  print  data  for  the  subsystems  (chunks  of  ele¬ 
ments)  specified  in  the  Preprocessor.  Will  not  compute  if  left 
blank. 

0  will  print  all  explosive  element  data  if  ECHECK  <  1000. 

1  will  print  only  those  explosive  elements  which  have  been  fully 
detonated  if  ECHECK  <  1000. 

2  will  not  print  any  explosive  element  data. 

Interval  of  cycles  at  which  cycle  data  will  be  printed.  If  NDATA 
2,  cycle  data  will  be  printed  for  every  other  cycle  (2,  4,  6,  etc.). 
If  left  blank,  cycle  data  will  be  printed  for  every  cycle. 

1  will  give  a  printer  plot  of  the  node  geometry  if  NABOR  nodes 
are  used. 
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INPUT  DATA  FOR  THE  POSTPROCESSOR 


The  function  of  the  Postprocessor  is  to  provide  plots  of  the  results  in  the  form  of  state 
plots  and  time  plots.  The  state  plots  show  results  for  the  entire  system  at  a  specified 
time,  and  the  time  plots  show  results  for  a  specified  variable  as  a  function  of  time. 

State  Plots 

Input  data  for  state  plots  are  summarized  in  Figure  6.  Included  is  the  capability  to 
plot  geometries,  velocity  vectors,  contours  of  several  variables,  various  jet  character¬ 
istics  and  behind  target  debris.  Plots  can  be  requested  in  the  order  of  increasing 
time  and  cycle  numbers.  The  plots  can  be  requested  by  either  time  or  cycle  number. 
By  using  the  time  option  it  is  possible  to  request  plots  without  having  access  to  the 
output  from  the  Main  Routine.  The  times  at  which  data  are  requested  must  simply 
coincide  with  those  specified  on  the  Data  Output  Cards  of  the  Main  Routine.  End 
state  plot  data  with  a  blank  card. 

Geometry  Plot  Card  (215,  F10.0,  412,  22X,  A30)  - 

1  Code  to  specify  geometry  plot. 


CYCLE 


Cycle  number  of  the  plot  which  is  desired.  The  cycle  numbers  of 
the  data  written  on  the  restart  tape  are  given  in  the  printed  out 
put  of  the  Preprocessor  and  the  Main  Routine.  If  CYCLE  -  0  the 
plots  are  requested  on  the  basis  of  time. 


TIME 


Time  of  the  plot  which  is  desired.  Plots  can  be  requested  bv 
either  TIME  or  CYCLE. 


0  will  use  the  r  and  z  axes  from  the  previous  plot.  This  option  al¬ 
lows  deformed  geometry  to  be  plotted  together  with  contours  or 
velocity  vectors,  for  instance. 


AXES 


1  will  automatically  compute  the  r  and  z  axes  to  include  all 
nodes.  The  vertical  axis  is  specified  to  be  10  units,  and  the' 
horizontal  axis  is  as  required,  using  the  same  scale  as  the  ver¬ 
tical  axis. 


ORIENT 


2  will  read  the  coordinate  limits  of  the  plot. 

S' 

0  will  specify  the  z  axis  as  the  vertical  axis  and  the  r  axis  as  the 
horizontal  axis. 

I  will  specify  the  r  axis  as  the  vertical  axis  and  the  z  axis  as  the 
horizontal  axis. 
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0  will  plot  the  grid  on  the  actual  r  coordinates. 

=<  1  will  plot  the  grid  on  the  negated  r  coordinates  only. 

2  will  plot  the  grid  on  both  the  actual  and  negated  r  coordinates. 
^3  will  not  plot  the  grid. 

0  plots  no  outline  around  the  edges. 

1  plots  an  outline  on  the  negated  r  coordinate  only. 


EDGE 


2  plots  an  outline  on  both  the  actual  and  negated  r  coordinates. 
-<  3  plots  an  outline  on  the  actual  r  coordinates  only. 


-1  is  same  as  EDGE  =  1  except  outlines  are  also  included  be 
tween  different  materials. 

-2  is  same  as  EDGE  =  2  except  outlines  are  also  included  be¬ 
tween  different  materials. 

-3  is  same  as  EDGE  =  3  except  outlines  are  also  included  be- 
,  tween  different  materials. 


=  1  will  plot  star  in  center  of  fractured  (partially  failed  element). 
Will  not  plot  star  if  left  blank. 

'O  will  not  plot  individual  node  points. 

1  will  plot  node  points  on  negated  r  coordinates  only. 


2  will  plot  node  points  on  both  the  actual  and  negated  r  coordi 
nates. 


NODE 


3  will  plot  node  points  on  actual  r  coordinates  only. 

=\  -1  is  same  as  NODE  =  1  except  projectile  nodes  are  drawn  as  a 
plus  sign  and  target  nodes  are  drawn  as  a  diamond. 


-2  is  same  as  NODE  =  2  except  nodes  are  drawn  as  symbols  in 
stead  of  points. 

-3  is  same  as  NODE  -  3  except  nodes  are  drawn  as  symbols  in 
stead  of  points. 


TITLE 


Title  printed  on  the  plot. 


Plot  Limits  Cards  for  Axes  =  2  (4F10.0)  —  This  card  specifies  the  portion  of  the 
problem  which  is  plotted.  Regions  beyond  those  specified  are  not  plotted.  The  ver 
tical  axis  is  scaled  to  10  units  and  the  horizontal  is  as  required.  The  scale  factor  used 
will  be  a  multiple  of  1,  2,  3,  5,  or  8  per  axis  unit. 

RMAX  =  Maximum  r  coordinate  included  in  the  plot. 

RMIN  =  Mimimum  r  coordinate  included  in  the  plot. 

/.MAX  =  Maximum  z  coordinate  included  in  the  plot. 

ZMIN  =  Minimum  z  coordinate  included  in  the  plot. 

Velocity  Vector  Plot  Card  (215,  F10.0,  512,  10X,  F10.0,  A30)  —  Only  previously 
undefined  variables  are  defined.  See  Geometry  Plot  Card  description  for  others. 

2  =  Code  to  specify  velocity  vector  plot. 

ARROW  =  1  will  place  arrowheads  on  the  velocity  vectors.  Leave  blank  for 

no  arrowheads. 

VSCALE  =  Velocity  which  will  give  a  velocity  vector  which  has  a  length  of 

1.0  using  the  scale  of  the  plot.  If  left  blank,  VSCALE  will  auto 
matically  be  determined  to  give  the  longest  vector  a  length  of 
two  percent  of  the  length  of  the  vertical  axis. 

Plot  Limits  Card  for  AXES  =  2  (4F10.0)  —  This  card  is  the  same  as  described  for 
the  geometry  plots. 

Contour  Plot  Card  (215,  F10.0,  712,  16X,  A30)  —  This  card  requests  contour  plots 
of  element  variables.  Contours  are  determined  by  first  computing  the  variable  quan 
titles  at  the  nodes  (i.e.,  the  nodal  pressure  is  the  average  of  the  pressures  of  all  ele 
ments  which  contain  the  node).  Then  the  contours  are  drawn  through  the  nodal 
quantities.  Only  previously  undefined  variables  are  defined.  See  Geometry  and  Ve¬ 
locity  Plot  Cards  for  others. 

TYPE  -  Code  to  specify  which  variable  is  requested.  Must  be  in  the 

range  of  3-10.  See  Figure  6  for  description  of  variables. 

NLINE  =  Number  of  contours  to  be  plotted.  If  NLINE  =  0,  six  contours 

will  be  plotted  at  values  of  5,  20,  40,  60,  80  and  95  percent  of  the 
range  between  the  minimum  and  maximum  variable  quantity 
limits. 
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=  Increment  at  which  symbols  are  placed  on  contour  lines.  SYM 
BOL  -  1  will  place  symbols  at  the  forward  end  of  each  contour 
line  within  an  element,  and  SYMBOL  =  5  will  place  symbols  at 
the  forward  end  of  every  fifth  element,  etc.  SYMBOL  =  0  will 
place  only  one  symbol  on  the  line. 

PRINT  =  1  will  print  the  nodal  quantities  of  the  specified  variable  on  the 

output  of  the  Postprocessor.  Leave  blank  for  no  printed  nodal 
output. 

Plot  limits  Card  for  AXES  =  2  (4F10.0)  —  This  card  is  the  same  as  described  for 
the  geometry  plots. 

Contour  Specification  Card(s)  (8F10.0)  —  Used  only  for  NLINE  >  0  on  Contour 
Plot  Card. 

Cl  ...  CN  =  Magnitude  of  contours  to  be  plotted 

One-Dimensional  Plot  Card  (215,  F10.0,  12,  10X,  12,  A30  —  For  one  dimensional 
geometry,  variables  are  plotted  as  a  function  of  the  z  axis.  The  plot  axes  are  divided 
into  ten  units  each.  Plot  types  must  be  in  the  range  of  2-10  as  shown  in  Figure  6.  All 
input  variables  have  been  previously  defined. 

Plot  Limits  Card  for  AXES  =  2  (4F10.0)  — 

VMAX  =  Maximum  value  of  the  dependent  variable  included  in  the  plot 

VMIN  =  Minimum  value  of  the  dependent  variable  included  in  the  plot 

ZMAX  =  Maximum  z  coordinate  included  in  the  plot 

ZMIN  =  Minimum  z  coordinate  included  in  the  plot 

Jet  Characteristics  Plot  Card  (215,  F10.0,  212,  26X,  A30)  —  This  card  requests 
plots  of  various  jet  characteristics  for  explosively  formed  penetrators,  hemi  charges 
and  shaped  charges.  Only  data  from  the  projectile  nodes  and  masses  are  plotted. 
This  allows  meaningful  plots  to  be  obtained  before  the  explosive  gases,  etc.,  are  dis 
carded.  These  plots  have  a  vertical  axis  of  10  units  and  a  horizontal  axis  of  10  units. 
Only  previously  undefined  variables  are  defined.  See  Geometry  Plot  Card  for  others. 

TYPE  Code  to  specify  which  plot  is  requested.  Must  be  in  the  range  of 

11-17.  See  Figure  6  for  description  of  the  options.  Plots  for  one¬ 
dimensional  runs  may  also  be  obtained  for  TYPE  -  11-15. 

0  will  accumulate  mass,  momentum  and  kinetic  energy  from  the 
slower  end  of  the  jet 

1  will  accumulate  from  the  faster  end  (tip)  of  the  jet 


At  CUM 


Plot  Limits  Card  for  AXES  =  2  and  TYPE  =  11-15  (4F10.0)  — 


VMAX  -  Maximum  value  of  the  dependent  variable 

VMIN  -  Minimum  value  of  the  dependent  variable 

// Z  MAX  =  Maximum  value  of  the  independent  position  or  velocity  variable 

ZZ  MIN  =  Minimum  value  of  the  independent  position  or  velocity  variable 

Plot  Limits  Card  for  AXES  =  2  and  TYPE  =  16-17  (4F10.0)  — 

UMAX  =  Maximum  r  coordinate  included  in  the  plot 

RMIN  —  Minimum  r  coordinate  included  in  the  plot 

R  Z  MAX  =  Maximum  velocity  included  in  the  plot 

R/Z  MIN  =  Minimum  velocity  included  in  the  plot 

Behind  Target  Debris  Plot  Card  (215,  F10.0,  412,  12X,  F10.0,  A30)  -  This  card 
requests  plots  of  behind  target  momentum,  kinetic  energy,  mass,  or  nodal  velocity. 
Only  velocities  in  the  R  Z  plane  are  included.  The  distribution  plotted  can  be  either 
that  present  at  a  flat  collector  plate  or  the  angular  distribution  at  infinity.  The  col 
lector  plate  model  plots  momentum/unit  area,  kinetic  energy/unit  area,  mass/unit 
area,  and  net  nodal  velocity  versus  plate  position.  The  angular  distribution  model 
plots  momentum/sterradian,  kinetic  energy/sterradian,  mass/sterradian,  and  net 
nodal  velocity  versus  angle  from  the  z  axis.  Only  previously  undefined  variables  are 
defined. 

TYPE  -  Code  to  specify  which  plot  is  requested.  Must  be  in  the  range 

18  21.  See  Figure  6  for  description  of  the  options. 


OPTION 


1 1  will  plot  results  for  flat  collector  plate 
I  2  will  plot  results  for  angular  distribution 


)  0  will  use  a  cutoff  velocity  of  one-tenth  the  most  negative  z  velocity. 
1 1  will  read  the  cutoff  velocity  from  the  Plot  Limits  Card 


NB1N  =■  Number  of  collector  bins  to  sum  variables.  Must  be  in  range  of 

1-99.  NBIN  =  0  will  give  90  angular  bins  (of  one  degree  each)  for 
angular  distribution  model. 

|  0  will  use  all  the  nodes 

NPART  !  1  will  use  projectile  nodes  only 

f  2  will  use  target  nodes  only 

Z  coordinate  of  the  flate  collector  plate  for  OPTION  =  1. 
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ZPLATE 


Plot  limits  Card  (5F10.0) 

VN1AX  =  Maximum  value  of  the  dependent  variable.  If  VMAX  =  VMIN 

0,  then  the  values  of  VMAX  and  VMIN  will  be  computed  from 
the  data  and  the  vertical  scale  will  be  automatically  determined 

VMIN  -  Minimum  value  of  the  dependent  variable.  Automatically  de 

termined  when  VMAX  =  VMIN  =  0. 

HMAX  =  Maximum  r  coordinate  of  collector  plate  for  OPTION  =  1,  or 

maximum  angle  (degrees)  for  OPTION  =  2 

HMIN  =  Minimum  r  coordinate  of  collector  plate  for  OPTION  =  1,  or 

minimum  angle  (degrees)  for  OPTION  -  2. 

V-CUTOFF  -  Cutoff  velocity  for  VCUT  =  1.  Nodes  with  a  net  velocity  less 
than  the  cutoff  velocity  will  not  be  included  in  the  data  for  the 
plots.  For  VCUT  =  0,  V-CUTOFF  will  be  assigned  a  cutoff  veloc¬ 
ity  of  one-tenth  the  most  negative  z  velocity. 

Kxtrapolated  Geometry  Plot  Card  (215,  F10.0,  612,  18X,  A30)  —  This  option  al 
lows  the  user  to  obtain  extrapolated  geometry  plots  at  times  much  greater  than  were 
computed.  Similar  options  as  for  Geometry  Plot  Card  (TYPE  =  1).  Only  previously 
undefined  variables  are  defined. 

T  EXTRAP  =■  Extrapolated  time  for  which  the  geometry  is  desired.  Nodal  posi 
tions  are  based  on  straight  line  extrapolation  using  positions 
and  velocities  from  the  specified  cycle  (or  time). 

Time  Plots 

Input  data  for  time  plots  are  summarized  in  Figure  7.  System  Plot  Cards  should  be 
input  first,  followed  bv  Individual  Node  and  Individual  Element  Plot  Cards.  The  var¬ 
iables  are  plotted  as  a  function  of  time.  Program  is  currently  dimensioned  such  that 
system  data  can  be  plotted  for  a  maximum  of  500  time  increments,  and  the  node  and 
element  data  for  a  maximum  of  1000  time  increments.  The  plot  axes  are  divided  into 
ten  units  each.  End  with  a  blank  card. 


System  Plot  Cards  (215,  5X,  F5.0,  4F10.0,  A18)  —  These  cards  request  plots  of  the 
system  variables.  Each  plot  contains  data  for  the  projectile,  the  target  and  the  total 
system  (projectile  plus  target).  These  data  must  have  been  previously  written  on  the 
plot  tape  by  setting  SYS  1  on  the  Plot  Card  in  the  Main  Routine. 

TYPE  Code  describing  the  type  of  plot.  See  Figure  7  for  description  of 

type.  Must  be  in  range  of  1  through  17. 


10  will  automatically  select  coordinates  to  include  max  and  min 
values  of  variable  for  total  duration  of  time 
1  will  read  the  coordinate  limits  of  the  plot 

SCALE  =  Factor  by  which  the  variables  are  multiplied  before  plotting 

TMAX  =  Maximum  time  included  on  horizontal  axis  if  AXES  =  1  (time) 

TM1N  =  Minimum  time  included  on  horizontal  axis  if  AXES  =  1  (time) 

V'MAX  =  Maximum  variable  included  in  vertical  axis  if  AXES  -  1 

VMIN  =  Minimum  variable  included  on  vertical  axis  if  AXES  1 

TITLE  --  Title  written  on  the  plot 

Individual  Node  Plot  Cards  (315,  F5.0,  4F10.0,  A18)  —  These  cards  request  plots 
of  nodal  variables.  These  data  must  have  been  previously  written  on  the  plot  tape  by 
specifying  the  requested  nodes  on  the  Designated  Nodes  Card  in  the  Main  Routine. 
Only  previously  undefined  variables  are  defined. 

TYPE  =  Code  describing  the  type  of  plot.  See  Figure  7  for  description  of 

types.  Must  be  in  the  range  of  18  through  30  Note  that  accelera 
tion  data  (TYPE  =  24-26)  may  be  incorrect  for  sliding  surface 
and  rigid  body  nodes. 

NODE  =  Specific  node  for  which  plot  data  are  requested. 

Individual  Element  Plot  Cards  (315,  F5.0,  4F10.0,  A18)  —  These  cards  request 
plots  of  element  variables.  These  data  must  have  been  previously  written  on  the  plot 
tape  by  specifying  the  requested  elements  on  the  designated  Elements  Cards  in  the 
Main  Routine.  Only  previously  undefined  variables  are  defined. 

TYPE  =  Code  describing  the  type  of  plot.  See  Figure  6  for  description  of 

types.  Must  be  in  the  range  of  31  through  39.  The  only  data 
available  for  explosive  elements  are  pressure  (TYPE  =  31)  and 
Internal  Energy  (TYPE  =  38). 

ELK  =  Specific  element  for  which  plot  data  are  requested 

PROGRAM  STRUCTURE  AND  FILE  DESIGNATION 

A  hierarchy  chart  for  the  Preprocessor  and  Main  Routine  is  shown  in  Figure  16.  The 
Postprocessor  for  state  plots  and  the  Postprocessor  for  time  plots  are  two  separate 
programs  The  corresponding  hierarchy  charts  are  shown  in  Figures  17  and  18.  The 
grouping  of  subroutines  is  given  in  Table  1  and  the  required  groups  at  subroutines 
Ibr  various  types  of  runs  are  given  in  Table  2 
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Figure  18.  Hierarchy  Chart  for  the  Time  Plots  Postprocessor 


The  file  designations  are  as  follows: 

ITAT  -  4  copy  of  Target  Addition  Cards  from  Input  file  when  Target  Ad 

dition  option  is  used. 

IN  -  5  Input  file 

IOUT  =  6  Output  file 


£ 


G 


IPLTIN  7  Restartable  plot  tape  for  time  plots,  which  is  read  by  the  Main 

Routine 

ITAPLT  -  8  Plot  tape  for  time  plots  written  by  the  Main  Routine 

ITAPIN  -  9  Restart  tape  ready  by  the  Main  Routine 

ITAPOT  -  9  Restart  tape  generated  by  the  Preprocessor  and/or  Main  Rou¬ 

tine.  When  ITAPOT  =  ITAPIN  the  restart  run  writes  data  on 
the  same  tape  from  which  it  read  the  restart  data.  Previous  cycle 
data  are  also  saved  on  the  tape.  If  ITAPOT  ^  ITAPIN,  the  re 
start  data  are  written  on  a  different  tape. 
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Table  1 .  Subroutine  Groupings  for  the  EPIC  2  Code 


EPIC2 

PREP 

MAIN 

NABOR 

P08T1 

P08T2 

EPIC2 

GEOM 

CHANQE 

CRUSH2 

POBT1 

POST2 

ECORE 

HEAT 

CHUNK 

INVAQ 

ACCUM 

EPLOT 

ECOUNT 

MATL 

CPCLCK 

MEGRU2 

AXIE8 

NPLOT 

E3HAPE 

MATLIB 

CRU8H 

NABOR2 

CONTUR 

PDATA 

EXTRA 

START 

DAMAGE 

SEARCH 

DAT8ET 

READTP 

HOATA 

ELOOP1 

SEEK2 

ECHOIN 

8YPLOT 

NBCORE 

ELOOP2 

STRE82 

EDGE 

NCORE 

EPIC1 

FORGE 

NFIX 

ERODE 

GEOM 

N3HAPE 

HEBURN 

MAP 

PPLOT 

LOAD 

MAPRZ 

PRINT 

LOOP 

ONQRAF 

RIQID 

MATCH1 

READCY 

SAVE 

MATCH3 

SIDE8 

SOATA 

MIEGRU 

TOP 

SORT 

NEWMAS 

VEL 

NLOOP 

V8 

PCALL 

PRES2 

PRES4 

QUAD 

RECALL 

REPORT 

SEEK1 

SLIOE 

STRE33 

TARGAD 

TPLOTS 

WINDOW 

INSTRUCTIONS  FOR  CHANGING  PROGRAM  DIMENSIONS 

The  dimensions  of  the  Preprocessor  and  Main  Routine  can  be  changed  by  re¬ 
dimensioning  the  arrays  in  common  blocks  NODE,  ELEMNT,  NABOR  and  MISC2 
The  corresponding  sizes  of  the  redimensioned  arrays  must  be  placed  in  the  DATA 
statement  in  the  calling  program,  EPIC-2.  A  description  of  the  array  sizes  for  the 
DATA  statement  is  given  in  subroutine  GEOM.  Some  additional  comments  are  as 
follows: 

•  If  no  elements  are  used,  the  ELEMNT  arrays  may  be  set  to  the  minimum  size 
MAXL  =  1. 

•  If  no  NABOR  nodes  are  used,  the  NABOR  arrays  may  be  set  to  the  minimum 
size  MAXNAB  =  1. 

•  If  NABOR  nodes  are  used,  the  NODE  array  sizes  and  NABOR  array  sizes 
should  be  identical  (MAXN  =  MAXNAB). 

CENTRAL  PROCESSOR  TIME  ESTIMATES 

For  various  problems  run  on  the  CDC  Cyber  176  computer  at  Eglin  AFB,  the  EPIC-2 
Main  Routine  has  used  0.0001  to  0.0002  central  processor  second  per  cycle  per  ele¬ 
ment.  The  explosive  elements  require  less  computing  time  than  the  solid  elements. 

EXAMPLE  PROBLEM 

Input  data  for  a  copper  projectile  impacting  a  steel  target  are  given  in  Figure  19. 
The  geometry  of  the  computed  response  is  shown  in  Figure  20  and  output  data  are 
shown  in  Figure  21. 


Figure  19.  Input  Data  for  the  Example  Problem 


SYSTEN  DATA  FOR  CYCLE  534  AT  TINE  *  0.400604E-04 


PROJECTILE 

TARGET 

TOTAL  SYSTEH 

NASS 

ERODED  NASS  (TOTALLY 

FAILED) 

0.131410E-02 

0.872645E-03 

0.206968E-01 

0.131526E-02 

0.220109E-01 

0.218791E-02 

TOTAL  ENER6Y 

KINETIC  ENERSY 
INTERNAL  ENER6Y 
PLASTIC  UORK 

0.349938E407 
0. 320827E+07 
0. 2911 11 E+06 
0.668562E+05 

0.247590E*07 
0.140229(407 
0.107361(407 
0. 833626E+06 

0.597528E407 

0.461056E407 

0.136472E407 

0.900483E406 

NAKINUN  R  COORDINATE 
NININUN  R  COORDINATE 
NAXIHUN  Z  COORDINATE 
NININUN  Z  COORDINATE 

0.  21 2466E+01 
O.OOOOOOETOO 
0.286858E401 
“0. 246614E+01 

0.302742E401 

0.926403E-01 

0.122743E401 

-0.251750E401 

0.302742E401 
0.000000 E400 
0.286858E401 
-0.251750E401 

C6  IN  R  DIRECTION 

C 6  IN  Z  DIRECTION 

0.000000£*00 

-0.193553E401 

O.OOOOOOETOO 

-0.567819E400 

O.OOOOOOETOO 

-0.649475E400 

LINEAR  NONENTUN  IN  R  DIRECTION 
LINEAR  NONENTUN  IN  Z  DIRECTION 
ANCULAR  NONENTUN  (CLOCKWISE  IS 

4)  « 

O.OOOOOOETOO 

-0.812704£402 

O.OOOOOOETOO 

O.OOOOOOE+OO 

-0.501396E402 

O.OOOOOOETOO 

0.000000E400 
-0. 1 31 410E403 
O.OOOOOOETOO 

LINEAR  VELOCITY  IN  R  DIRECTION 
LINEAR  VELOCITY  IN  Z  DIRECTION 
NET  LINEAR  VELOCITY 

ANSULAR  VELOCITY  (CLOCKWISE  IS 

♦  >  « 

0. OOOOOOE+OO 
”0.  61 84498*05 
0.61 8449ET05 
O.OOOOOOETOO 

O.OOOOOOETOO 
-0. 242257(404 
0.2422S7E404 
0.000000E400 

O.OOOOOOETOO 

-0.397022E404 

0.597022E404 

O.OOOOOOETOO 

CYCLE  534  DATA  SAVED  ON  FILE  9  FOR  RESTART  AND/OR  PLOTS 


Figure  21 .  Output  Data  for  the  Example  Problem 


SECTION  III 

CONCLUSIONS  AND  RECOMMENDATIONS 


User  instructions  have  been  provided  for  the  EPIC-2  code.  This  code  can  be  used  foi 
a  wide  range  of  problems  involving  high-velocity  impact,  explosive-metal  inter 
action,  and  analysis  of  various  material  characterization  tests. 
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